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ABSTRACT 


The  general  problem  of  the  response  of  a  cantilevered  beam  to  flow  over  its 
surface  is  considered  experimentally  and  theoretically.  The  measured  flow-induced 
modal  vibratory  motion  of  a  nonsinging  beam  is  compared  to  theoretical  estimates 
of  inflow  turbulence  excitation  and  boundary-layer  excitation.  The  comparison  in¬ 
dicates  that  while  the  response  to  turbulent  inflow  is  dominant  at  low  frequencies, 
the  response  of  the  strut  to  its  own  boundary  layer  is  important  at  high  frequen¬ 
cies.  The  magnitude  of  hydrodynamically  induced  damping  is  also  characterized  ex¬ 
perimentally.  It  is  shown  that  results  agree  favorably  with  an  approximate  expres¬ 
sion  based  on  finite-aspect-ratio,  unsteady  airfoil  theory.  Loss  factors,  based  on 
entrained  mass,  are  found  to  be  inversely  proportional  to  a  reduced  frequency  based 
on  the  width  of  the  strut  and  inflow  speed.  Finally,  a  wind  tunnel  study  of  the 
statistical  properties  of  the  boundary  layer  formed  on  the  strut  is  described.  The 
results  disclose  that  flow  separation  at  the  leading  edge,  which  is  sensitive  to  angles 
of  attack,  generates  a  low-frequency  pressure  field  that  is  markedly  higher  than 
that  normally  encountered  in  boundary  layers.  At  high  frequencies,  the  pressure 
field  is  influenced  by  the  local  flow  parameters  normally  used  in  boundary-layer 
studies. 


ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  was  funded  by  the  Naval  Ship  Systems  Command  under  Sub- 
project  S-F35  452  007,  Task  15129,  Work  Unit  1942-057. 

1.  INTRODUCTION 

The  motion  of  a  beam  in  water  is  influenced  by  water  loading,  both  hydrodynamically 
and  acoustically.  In  the  case  of  an  unbaffled  beam  without  a  mean-fluid  motion,  the  water- 
loading  is  due  to  entrained-fluid  inertia  arising  from  fluid  forces  in  phase  opposition  with  ac¬ 
celeration,  to  acoustic  radiation  arising  from  compressible  fluid  forces  in  phase  opposition 
with  the  beam  velocity,  and  to  viscous  forces  arising  from  shearing  in  the  fluid  associated 
with  flow  around  edges  and  near  nodal  points.  The  viscous  effect  is  generated  as  the  beam 
moves  normal  to  its  neutral  surface  and  sets  up  a  local  eddy  field  which  is  dissipated  by  vis¬ 
cosity.  The  viscous  losses  have  been  studied  on  free-free  beams  by  Blake  (1972),1  who  has 
shown  the  dependence  of  this  damping  on  fluid  properties  and  beam  geometry.  Acoustic 


1  Blake,  W.  K.,  “On  the  Damping  of  Transverse  Motion  of  Free-Free  Beams  in  Dense  Stagnant  Fluids.”  Shock  and  Vibra¬ 
tion  Bulletin  Vol.  42,  Part  4,  pp.  41-55  (Jan  1972).  A  complete  listing  of  references  is  given  on  pages  66  and  67. 
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and  inertial  fluid  reactions  on  free-free  beams  have  been  studied  by  Blake  ( 1 97 1  )2  for  cases 
in  which  the  length  of  the  acoustic  wave  is  substantially  larger  than  the  beam  width,  thick¬ 
ness,  and  structural  wave  length.  Also,  for  beams  at  frequencies  lower  than  acoustic  coinci¬ 
dence,  the  experiments  of  Reference  2  have  shown  that  modal  coupling  may  be  easily 
neglected. 

When  there  is  mean-fluid  motion  parallel  to  the  plane  of  the  neutral  surface  of  the  beam, 
the  hydrodynamic  forces  are  generated  by  the  potential  flow  around  the  beam.  Superposition 
of  the  steady  inflow  with  the  vibration  velocity  of  the  beam  is  manifested  in  an  oscillating 
angle  of  attack  which  results  in  fluid  forces  and  moments.  The  motion  of  the  beam  gener¬ 
ates  vorticity  in  its  wake,  and  the  fluid-reaction  force  is  equal  to  the  time  rate  of  change  of 
momentum  imparted  to  the  wake.  The  relative  importance  of  the  forces  and  moments  in¬ 
volved  in  the  circulation  fluctuations  depend  on  the  mode  of  vibration;  they  are  also  capable 
of  introducing  coupling  between  the  different  modes.  We  will  be  concerned  with  cases  in¬ 
volving  relative  magnitudes  of  entrained  mass,  which  are  less  than  or  comparable  to  the  mod¬ 
al  mass  of  the  strut.  Also,  our  considerations  will  be  restricted  to  low  enough  speeds  that 
the  structure-fluid  instabilities  leading  to  flutter  divergence  are  not  important.  At  these  low 
speeds  and  low  ratios  of  added  mass-to-structure  mass,  the  hydrodynamic  reactions  are  large¬ 
ly  restricted  to  forces  and  moments  which  oppose  the  beam  motion  and  are  proportional  to 
the  beam  velocity;  they  increase  the  observed  damping  of  the  strut. 

In  addition  to  the  reaction  fluid  forces,  we  must  characterize  the  flow-induced  excita¬ 
tion  forces.  These  can  be  caused  by  the  boundary  layer  formed  on  the  strut,  by  the  forma¬ 
tion  of  the  viscous  wake  at  the  trailing  edge,  and  by  interaction  with  turbulent  velocity  fluc¬ 
tuations  in  the  inflow.  In  the  case  of  boundary -layer  excitation,  the  local  pressure  field  is 
more  or  less  distributed  over  the  width  and  span  of  the  beam.  This  excitation  will  be  altered 
by  changes  in  the  angle  of  attack  of  the  beam  and,  possibly,  by  beam  vibration  if  the  modal 
amplitudes  are  high  enough. 

When  wake-vortex  excitation  occurs,  the  pressures  are  localized  at  the  trailing  edge,  and 
they  are  often  characterized  by  discrete-frequency  excitation.  The  vortex-induced  excitation 
is  most  intense  when  the  separation  points  of  the  boundary  layer  are  well  defined. 

A  third  source  of  excitation  involves  the  reaction  of  the  beam  to  turbulence  in  the  in¬ 
flow.  These  pressures  are  induced  because  of  the  unsteady  local  angle  of  attack  generated 
by  the  inflow  turbulence.  The  pressures  are  maximum  at  the  leading  edge.  A  recent  inves¬ 
tigation  by  Mugridge  (1970), 3  which  is  an  extension  of  the  early  theoretical  work  of 


^Blake,  W.  K.,  “Radiation  from  Free-Free  Beams  under  Influences  of  Light  and  of  Heavy  Fluid-Loading,”  NSRDC  Report 
3716  (Nov  1971). 

^Mugridge,  B.  D.,  “The  Generation  and  Radiation  of  Acoustic  Energy  by  the  Blades  of  a  Sub-Sonic  Axial  Flow  Fan  Due 
to  Unsteady  Flow  Interaction,”  University  of  Southampton  (England)  Ph.  D.  Thesis  (Jan  1970). 


Reissner  (1947), 4  characterizes  the  unsteady  forces  affecting  an  airfoil  in  terms  of  the  spec¬ 
trum  of  inflow  turbulence.  The  results  are  expected  to  apply  in  the  current  investigation, 
except  when  effects  of  finite  thickness  become  important. 

The  purpose  of  the  investigation  was  to  measure  the  flow-induced  vibration  and  damp¬ 
ing  of  a  cantilevered  strut  in  water  and  to  compare  those  findings  to  theoretical  predictions. 
To  assist  in  the  theoretical  predictions,  a  subsidiary  wind  tunnel  experiment  was  performed 
to  characterize  the  statistics  of  the  boundary-layer  flow  field,  developed  on  a  rigid,  two- 
dimensional  strut  of  geometrically  similar  cross  section.  The  air  experiment  was  performed 
at  Reynolds  numbers  based  on  a  chord  near  to  those  experienced  in  the  water  test.  This 
measurement  provided  necessary  scaling  criteria  for  use  in  estimates  of  the  response  of  a  flex¬ 
ible  strut  to  flow.  Also  in  this  report,  estimates  of  the  dynamical  responses  of  the  beams  to 
inflow  are  compared  with  measurements.  A  theoretical  discussion  precedes  a  description  of 
the  measurements  in  order  to  establish  normalization  factors  with  which  to  examine  experi¬ 
mental  results. 


2.  AN  ANALYTICAL  DESCRIPTION  OF  THE  MOTION  OF  A  BEAM  IN  WATER 
FORMULATION  OF  THE  BEAM-VIBRATION  FIELD 

We  will  consider  the  response  of  a  beam  to  random  hydrodynamic  pressure  fields  which 
are  generated  on  both  sides  of  the  beam  by  any  of  the  three  sources  outlined  previously.  The 
fluid  loading  in  the  forms  of  acoustic,  inertial,  and  unsteady  lift  pressures  are  included  in  the 
analysis.  The  basic  analysis  essentially  follows  the  format  of  the  Lin  (1967)5  normal  mode 
analysis  of  a  linear  system  when  the  statistical  properties  of  the  excitation  are  specified. 
Leehey  (1968)6  and  Davies  (1969)7  have  used  the  same  approach  to  estimate  the  response  of 
panels  to  boundary-layer  turbulence.  The  relationships  which  are  presented  here  could  be 
straightforwardly  derived  by  modifying  any  of  the  results  available  in  the  previously  described 
references.  For  completeness,  however,  the  basic  development  has  been  outlined  as  follows. 
The  beam  motion  is  a  low-frequency  limit  for  the  motion  of  a  cantilevered  plate  for  which 
three  sides  are  free,  and  the  fourth  side  is  clamped.  In  the  experiments  described  later,  both 
beam  and  plate  modes  were  encountered. 


^Reissner,  E.,  “On  the  General  Theory  of  Thin  Airfoils  for  Non-Uniform  Motion,”  National  Advisory  Committee  for 
Aeronautics  TN  946  (Aug  1947). 

^ Lin,  Y.K.,  “Probablistic  Theory  of  Structural  Dynamics,”  McGraw-Hill,  Inc.,  New  York  (1967). 

^Leehey,  P.,  “Trends  in  Boundary  Layer  Noise  Research,”  Proceedings  of  the  Air  Force  Office  of  Scientific  Research- 
University  of  Toronto,  Institute  for  Aerospace  Studies  Symposium,  Toronto,  Ontario,  Canada  (20-21  May  1968). 

^Davies,  H.  G.,  “Sound  from  Turbulent  Boundary  Layer  Excited  Panels,”  Massachusetts  Institute  of  Technology  Acoustics 
and  Vibration  Laboratory  Report  70208-2  (Feb  1969). 
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Figure  1  shows  the  strut  under  consideration;  it  is  positioned  so  that  the  flow  is  paral¬ 
lel  to  the  neutral  surface  and  perpendicular  to  the  longest  dimension  L.  The  beam  has  width 
w  and  thickness  h.  Space  coordinates  in  the  plane  of  the  neutral  surface  are  designated  by 
xT;  time  is  designated  by  t.  Fluctuating  hydrodynamic  pressures  pf  (x",t)  are  generated  on 
each  surface,  and  the  pressures  on  the  upper  and  lower  surfaces  may  or  may  not  be  corre¬ 
lated  with  each  other,  depending  on  the  manner  of  flow  excitation.  Fluid-reaction  pressures 
are  generated  by  the  blade  motion  so  that  those  on  the  upper  and  lower  surfaces  are  corre¬ 
lated.  In  these  cases,  we  let  ph(?T,  t)  denote  the  differential  of  the  flow-dependent  hydro- 
dynamic  reaction  pressures  associated  with  fluid  lifts  and  moments  on  the  surfaces,  and  we 
let  pa(x',t)  denote  the  differential  of  flow-independent  acoustic  and  inertial  reaction  pres¬ 
sures  on  the  upper  and  lower  surfaces.  We  will  neglect  viscous  damping  in  the  theoretical 
analysis.  The  equation  of  motion  for  the  vibration  velocity  v(x,  t)  of  the  plate  of  uniform 
thickness  and  density  is 

DV4  v(x\t)  +  pp  h(3  v  (x",t)  +  pp  h  v(x,t)  =  pf  (x’.t)  -  ph  (x'.t)  -  pa  (x’.t)  (1) 

where  pp  is  the  volume  density  of  the  strut  material,  and  the  dots  above  symbols  denote  time 
derivatives.  The  stiffness  is 


!  Eh3 
D  = - 

12  1  -  (j2 

where  E  is  Young’s  modulus,  and  o  is  Poisson’s  ratio  of  the  strut  material.  Mechanical  damp¬ 
ing  caused  by  friction  in  the  clamp,  and  hysteretic  material  losses  are  represented  by  an  ad 
hoc  mechanical  damping  coefficient  (5 .  The  strut  velocity  can  be  decomposed  into  its  gener¬ 
alized  time  Fourier  components 


V(x>)  =  — 

l  7T 


which  are  expandable  in  the  normal  modes  of  the  strut  ^mn(x  ) 


V(x  ,co)  =^Vmn(w)  *mn  (x) 

m,n 

where 


fw/2  r  L  _ 

J  J  (x)dx  =  L  •  w 

-w/2  J0 


oo 

j"  V()t,t)  eiwt  dt 

—  oo 


(2) 


(3) 


4' 


In  general,  ^mn  (x )  can  be  separated  into  eigenfunctions  which  are  dependent  on  the  coor¬ 
dinates  x  and  z  in  Figure  1 ,  i.e. 


(4) 


where  \pm  (z)  and  0n(x)  are  separate  eigenfunctions  describing  m-th  order  axial  variation  and 
n-th  order  lateral  variation  of  the  (m,n)  mode  pattern.  For  slender  beams  we  can  simplify 
our  discussions  by  considering  separately  those  modes  which  possess  node  lines  either  paral¬ 
lel  or  perpendicular  to  the  long  axis  of  the  beam.  The  former  family  of  modes  will  be  called 
torsion  modes  or  (o,n)  modes;  the  latter  family  will  be  called  bending  modes  or  (m,o)  modes. 
The  modes  that  have  a  combination  of  these  nodal  patterns,  i.e.,  (m,n)  modes,  will  be  called 
plate  modes.  The  plate  modes  were  not  often  encountered  in  the  experimental  work,  which 
will  be  discussed  later. 

Mode  shapes  for  modes  of  cantilevered  plates  have  been  theoretically  calculated  as  well 
as  measured;  a  review  of  the  results  of  many  investigations  has  been  given  by  Leissa  (1969).8 
Letting  m  equal  the  number  of  node  lines  along  the  width  and  n  be  the  number  of  node  lines 
along  the  length,  the  modes  of  a  beam  can  be  approximated  by  the  following  closed  forms* 


*0,1  &  =  ”  v'5"  [w) *  1  (kiz> 

*0,2  <X  >  =  <Mk2X)  [>Mk0Z)  +  ^  MklZ> 


(5A) 


(5B) 


*i,l  <x)  = 


2x 


w 


1  2 
-  (k„z)  +  i (kjZ)  +  -  ij/2  (k2z) 


(5C) 


*m,o(X)=  (kmZ>  (5D) 

where  i pm  (km  z)  is  the  orthogonal  set  of  mode  shapes  for  a  cantilevered  beam  with  m  node 
lines  normal  to  its  length,  and  <pn  (knx)  is  the  mode  shape  for  a  free-free  beam  with  n  node 
lines  normal  to  its  width.  The  modal  wave  numbers  are  given  by 


8Leissa,  A.  Q.,  “Vibration  of  Plates,”  National  Aeronautics  and  Space  Administration  Report  SP-160  (1969). 

♦Inspection  will  show  that  Equations  (5)  are  only  approximately  normalized  according  to  Equation  (3).  This  is  due  to 
the  series  truncations  necessary  to  obtain  the  simplified  forms.  Lack  of  normalization  is  most  important  for  the  (1,1)  mode. 
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km  =  (2m  +  1)  tt/2L  m  >  0 

and 

kn=0  0  <  n  <  2 

kn  =  (2n  -  1)  7t/2w  n  >  2 

The  generalized  time  Fourier  transform  of  Equation  (1),  with  subsequent  substitution  of 
Equation  (2)  into  Equation  (1)  gives,  using  Equation  (3) 


[° kmn  +iMpph0-  Pphco2  ]  Vmn(w) 


ico 

Lw 


fw/2  r  L 

J  J  [Pf 

_  /-I  •'A 


-w/2  •'O 


(x  ,co)  -Pa  (x  ,co)  -  Ph  (x  ,co)]  *mn  (x  )  dx 


(6) 


where  k2  =  k2  +  k? 

mn  m  n 

If  we  make  the  advantageous  assumption  of  acoustically  uncoupled  modes,  we  can  write 
the  acoustic  back  reaction  as  a  summation  over  modes,  i.e. 

Pa  (x ,  co)  -  ^7(rmn  +  i«  mmn )  >Tmn  (x  )  Vmn  (w)  (7) 

m,n 


where,  for  the  (m,n)  mode,  rmn  is  the  radiation  resistance  per  unit  area,  and  mmn  is  the 
entrained  fluid  mass  per  unit  area.  For  free-free  beams,  Blake2  has  shown  that 


for  ka/km  «  1,  and  kaw  <  4,  where  ka  is  the  acoustic  wave  number.  For  a  clamped-free 
beam,  Equation  (8)  is  likely  to  be  of  the  correct  decimal  order  of  magnitude.  The  entrained 
inertia  per  unit  area  for  beams  has  been  given  by  Blake  for  ka/km  «  1  to  be 


mm,0=Po7W  kmW<  1  (9A> 

mm,0  =P0  f  +kaw)_1  kmw>1  (9B) 
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Both  of  these  relationships  are  expected  to  apply  approximately  for  all  of  the  (m,0)  modes 
considered  in  the  experiments.  For  modes  involving  torsion,  rmn  is  likely  to  be  less  than  that 
predicted  using  Equation  (7).  In  using  Equation  (7),  we  have  assumed  that  the  fluid  imped¬ 
ance  is  locally  proportional  to  the  transverse  velocity  of  the  strut.  In  the  special  case  of  the 
torsion  (0,1)  mode,  the  fluid  inertial  impedance  is  manifested  in  an  entrained  moment  of 
inertia,  which  can  be  approximated  using  two-dimensional  potential  flow  relationships.  This 
special  case  will  be  further  examined  at  the  end  of  this  section. 


HYDRODYNAMIC  COUPLING  AND 
DAMPING  EFFECTS 

When  there  is  no  mean  flow  past  the  beam,  the  important  fluid  reactions  to  motion  are 
acoustic  and  inertial.  The  imposition  of  a  mean  flow  parallel  to  the  plane  of  the  neutral  sur¬ 
face  (Figure  1 )  causes  additional  forces.  Given  the  Fourier  component  of  vertical  velocity  on 
the  beam  surface  V(x,  to),  the  local  small  angle  of  attack  to  the  flow  of  mean  velocity  UQ  is 
V (x,w)/U0 .  This  Fourier  component  of  angle  of  attack  generates  a  generalized  unsteady 
pressure  having  a  major  component  in  phase  opposition  with  V(x  ,co).  Depending  on  whether 
the  motion  is  characterized  by  bending  (m,0)  modes  or  by  torsion  (0,n)  modes,  this  unsteady 
differential  pressure  field  ph  0?,co)  has  a  chordwise  distribution  which  can  provide  different 
relationships  between  net  forces  and  moments. 

The  basic  two-dimensional  theory  for  the  oscillating  hydrofoil  has  been  worked  out  in¬ 
dependently  by  Theodorsen  ( 1 9 3 5 )9  and  by  von  Karman  and  Sears  (1938).10  The  theoreti¬ 
cal  expressions  for  the  more  appropriate  three-dimensional  theory  have  been  derived  in  closed 
form  by  Reissner.4  The  expressions  are  unwieldly  for  our  purposes,  so  we  will  make  use  of 
certain  closed-form  approximations  to  numerical  solutions  of  the  Reissner  expressions.  These 
solutions  were  obtained  by  Lawrence  and  Gerber  (1952)11  and  were  experimentally  verified 
by  Laidlaw  and  Halfman  (1956).12  For  the  beam  (m,0)  modes  of  our  problem,  a  three- 
dimensional  correction  to  the  two-dimensional  expressions  will  be  a  function  of  reduced  fre¬ 
quency  gjc/U0,  where  c  is  the  half-width  of  the  beam,  and  local  aspect  ratio  =  7r/wkm. 
Since  the  lift  and  moment  distributions  are  not  uniform  over  the  chord  of  the  strut,  the  flow 
can  induce  modal  coupling  as  indicated  in  the  following  analysis. 


^Theodorsen,  T.,  “General  Theory  of  Aerodynamic  Instability  and  the  Mechanism  of  Flutter,”  National  Advisory  Com¬ 
mittee  for  Aeronautics  Report  496  (1935). 

l^von  Karman,  T.  and  W.  R.  Sears,  Journal  of  the  Aeronautical  Sciences,  Vol.  5,  p.  379  (1938). 

1 1  Lawrence,  H.  R.  and  E.  H.  Gerber,  Journal  of  the  Aeronautical  Sciences,  Vol.  19,  p.  769  (1952). 
l^Laidlaw,  W  R.  and  R.  L.  Halfman,  Journal  of  the  Aeronautical  Sciences,  Vol.  23,  p.  117  (1956). 
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The  local  hydrodynamic  pressure  ph  (x  ,<o)  can  be  written 

p,  (tw)  =  2  <n>  •  2«o  (10) 

m  ,n 

where  we  have  assumed  that  the  pressure  is  locally  proportional  to  the  vertical  velocity;  how¬ 
ever,  it  is  weighted  by  a  chordwise  lift-distribution  function  £mn(?)  that  depends  on  the  mode 
of  motion.  The  fluid  density  is  p0  ;  Cmn  (12),  where  12  =  coc/UQ ,  we  will  call  the  modal 
oscillatory  lift  coefficient,  which  is  a  function  of  mode  type  and  of  aspect  ratio  7r(k  w)-1. 
Equation  (10)  can  be  given  more  commonplace  significance  if  we  consider  the  special  case  of 
two-dimensional  translatory  motion.  This  case  is  approximated  in  practice  by  the  limit 
km  =  k0  -*■  0.  The  total  lift  distribution  dL/dz  for  a  two-dimensional  strip  is  given  by  ' 


dL 

dz 


f  w/2 
-w/2 


ph  (x,z,co)dx  =  2ir  pQ  wUQ  C1>0  (12) 


f  1/2 

-1/2 


£io(x)d 


w 


(ID 


where  the  term  in  square  brackets  is  the  two-dimensional  Theodorsen  function  C  (£2). 9 
Alternatively,  for  the  (0,1)  mode  of  an  infinite  beam,  we  have  two-dimensional  pitch  motion 
for  which  the  spanwise  moment  distribution  is 


dM 

dz 


f  w/2 

7T  ,  T  2v/3~ 

~  f1'2  X  2  x" 

-w/2 

Ph  (x>z>w)  xdx  =  fg  ^0  w  uo  — 

C01(12)  32£01(x)-  d- 

u’  i1/2  u’  W  W 

(12) 


where  the  term  in  square  brackets  is  [  1  -  C(12)  ]  ;  see  Sears  (1940).13  The  high-frequency 
limit  of  C(12)  is  Lim  C(12)  =  1/2. 

If  we  substitute  Equations  (7)  and  (10)  into  Equation  (6),  we  can  see  that  the  effect  of 
£mn  (x)  Is  t°  couple  modes,  if  £mn  (x)  is  not  independent  of  x.  We  obtain 


z“  (co)  Vmn  (co)  =  P  2 7T  C..  (12) pn  U„  V..  (co) 

mn  v  '  mn  v  '  mn  ijv  0  ij 

f  >.j 


_1_ 

Lw 


f  w/2  r  L 

J  J 

Z..r/n  •'a 


-w/2  "0 


(x)  ^y(x,z)^mn(x,z)dxdz 


(13) 


13 Sears,  W.  R.,  Journal  of  the  Aeronautical  Sciences,  Vol.  8,  p.  104  (1940). 


where 


The  (m,0)  modes  can  be  coupled  with  types  (m,l)  modes  because  a  heaving  motion  of  the 
beam  produces  a  moment  as  well  as  a  lift,  and  a  pitching  motion  of  the  beam  produces  a  lift 
as  well  as  a  moment.  The  first  two  factors  of  the  summation  term  of  Equation  (13)  can  be 
written  as 


(  ~  1  fw/2 

Pm ,0  =  i^P0  27r  J  Vm,0  (">  Cm,0  &)—  J  t  (x)  <A0  (x)  0Q  (x)  d  X 

h  l  W  -w/2 

1  fW/  2  ) 

+  Vm,l  Cm,l  J  £m,i  (X)01  (X)0  (X)dx  (14) 

W  -w/2  J 

where  we  have  used  Equation  (4).  To  evaluate  approximately  the  relative  sizes  of  the  modal 
amplitudes  Vm  0  (to)  and  Vml  (to),  we  note  that  Equation  (13)  can  be  evaluated  at  the  res¬ 
onance  frequency  of  the  (m,0)  mode  tom  0  for  U0  =  0  so  that 


Vm,0 


Here  we  have  assumed  that  the  Pmn  (to)  are  of  comparable  magnitude,  and  we  have  let 

f 


I  co  | 


^Ph  +  rm,0 

p„  h  +  mm  n 

rp  m  ,0 


where  r?  « 1  is  the  loss  factor  of  the  (m,0)  mode  in  still  water.  With  flow,  the  loss  factor 
increases,  but  if  it  remains  very  much  less  than  unity,  the  ratio  of  the  torsional  modal  ampli¬ 
tude  to  the  bending  modal  amplitude  is  still  likely  to  be  very  small  when  co  ,  and  to  „ 

m ,1  m  ,U 

are  sufficiently  different.  Furthermore,  in  the  case  of  £2  >  1 ,  C  (£2)  and  [  1  -  C  (£2)]  are 
both  approximately  one-half  so  that  both  coefficients  of  the  modal  amplitudes  in  Equation 
(14)  are  of  a  similar  decimal  order  of  magnitude.  Thus  for  light  damping,  for  which  r?  <  1, 
for  speeds  low  enough  that  £2  >  1,  and  for  com  j  considerably  different  than  tom  0  ,  this  type 
of  modal  coupling  can  be  neglected. 
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Equation  (1 3)  as  it  applies  to  beam  modes  now  becomes 


Dkm,0  -  w2(Pph  +  mm,o)  +  iM  trm,0+7rPoU0Cm,0  ^  Vm,o(W) 


=  i  OJ  Pm  Q  (oj)  (17) 

f’ 

The  lift  function  Cm  Q  (£2)  now  incorporates  the  integrated  lift-distribution  function;  it  still 
includes  a  three-dimensional  effect,  and  it  parallels  the  expression  of  Sears.13  The  lift  func¬ 
tion  accounts  for  spanwise  load  variation  but  neglects  the  overall  finite  length  L  of  the  beam. 
For  example,  for  large  L,  CQ  0  (£2)  is  just  the  Theodorsen  function.  Resonance  conditions  for 
the  lightly  damped,  clamped-free  beam  in  flow  occur  at  all  frequencies  for  which  Dk*  - 
(pp  h  +  mm )  =  0.  The  second  index  in  the  subscript  will  hereafter  be  dropped  when  beam 
modes  are  considered.  The  total  loss  factor  rjt  under  flow  conditions  is  then 


_  =  - - 

w(pp  h  +  rnm  ) 


PphcjT?m  +rm  +7rP0U0Cm,0<n) 


CO  (pp  h  +  mm  ) 


where  mm  is  given  by  Equation  (9), 

)3  has  been  replaced  by  coT?m  ,  and 
r?  is  the  mechanical  loss  factor. 

■m 

Equation  (18)  shows  that  as  UQ  is  increased,  the  hydrodynamic  damping  increases.  It  is 
pertinent  to  note  here  that  flow  effects  on  the  added  mass  are  neglected,  since  they  occur 
only  at  reduced  frequencies  less  than  unity.  Furthermore,  the  calculations  of  Lawrence  and 
Gerber1 1  show  that  the  three-dimensional  effects  of  structural  vibration  make  the  flow  con¬ 
tribution  to  entrained  inertia  even  more  insignificant. 

A  hydrodynamic  loss  factor  can  be  defined  in  terms  of  the  hydrodynamic  resistance 
and  entrained  inertia 


%  = 


7rp0  U0  Cm>0  (£2) 


For  km  w  <  1  the  entrained  mass  is  given  by  Equation  (9 A),  and  Cm  0  can  be  approximated 
by  its  two-dimensional  value  Cm  0  (£2)  —  1/2  for  £2  >  1,  and  7r/km  w  »  1,  so  that 


%  =  l/£2  (20) 

is  a  limiting  value  for  low-order  modes.  At  low-aspect  ratio,  i.e.,  7r/kn  w  <  1,  we  can  use  the 
Lawrence  and  Gerber11  small-aspect-ratio  approximation  to  give 
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(21) 


Cm,0  -  ff/2km  w 

and  combine  this  with  Equation  (9B)  to  give 


tjh  =  1/2  ft 

for  high-order  modes  at  high  reduced  frequencies. 

In  the  case  of  torsion  (0,n)  modes,  a  similarly  defined  hydrodynamic  loss  factor  is  ob¬ 
tained  by  using  the  Sears13  two-dimensional  moment  function.  The  flow-induced  resistance 
to  the  angular  velocity  of  the  section,  e.g.,  rQ  j ,  must  be  normalized  by  the  entrained  moment 
of  inertia  per  unit  area,  i.e.,  we  can  define  a  loss  factor  using  Sears13  expressions 


_  ro,i  ff/16p02  c2U0 

05  If  Joj  7r/128p0w3J  ^ 


(22) 


This  is  just  double  the  two-dimensional  loss  factor  induced  by  the  bending  motion.  The  ratio 
of  the  resonance  frequency  of  the  (0,1)  mode  in  water  to  that  in  air  is  just  the  square  root  of 
of  the  ratio  of  IB  to  Ifi  +  If,  where  IB  =  pph  •  c2  / 1 2.  Thus 


03 


wet 


CO. 


dry 


—i  1/2 


h  +  ^ 


(23) 


RESPONSE  OF  THE  WATERLOADED  BEAM 
TO  A  RANDOM  PRESSURE  FIELD 

We  will  assume  that  the  hydrodynamic  excitation  pressure  fields  on  each  side  of  the 
strut  are  uncorrelated  with  each  other.  Thus,  the  response  of  the  strut  to  large-scale  inflow 
turbulence  is  not  identically  given  by  the  following  analysis;  that  case  will  be  discussed  in 
Chapter  5.  The  generalized  time  transform  of  pf  (?,t)  is  a  stochastic  variable  as  well  as  its 
generalized  modal  coefficient  pm  0(co),  which  is  used  in  Equation  (17)  and  is  defined  follow 
ing  Equation  (13).  The  modal  velocity  coefficient  Vmn  (a>)  is  also  a  stochastic  variable.  We 
will  assume  statistical  stationarity  and  ergodicity  in  regard  to  time-averaging.  The  auto- 
spectral  density  of  the  vibration  velocity  measured  at  a  location  xQ  on  the  strut  is,  using 
Equation  (2), 

<t>w(oj,xQ )  =  2tt  <V(x0,w)  V*(x0A>)> 

=  2,  £ < <“>  Yt.  <">  >  *L  <70  ) 

m,n 

where  the  brackets  (<  > )  denote  an  ensemble  average. 
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Using  Equation  (17)  we  have  for  the  local  autospectral  density  of  velocity 


4>+}'X CO,X0)  =  2 


z 


V  n<*0>  I  f 


H- 


„  2  2  j  j  0p(X,y,Oj)'I'mn(x)'I'mn(y)dx  dV 

"  lzmnMI2  A?  ^A,  "A 


mnv  '  mn 

“s  *  ~s 

>iv  ;  J  V1'  '  Wi’  ')'$  y;.'L' 


i>|  .*7, *7,"  •'»yrrv  ».f*  u‘,*  iOY  i  ‘  ;  ■■■  -  ’>1":  ■  *  -'i  if r,  •  .^'.'.5  Ky  **i 

where  the  total  impedance  , 

ics/wi-no  'nbmi  rt- ■  i  ;  r;i  to_;Uix»vc  icfij&v.  v.:  • 

Zm.n;  ^),=  ::[D.kn,n  '  >  t.  >&]:  (iCj)"  ,•  W  :,iri,oi 

includes  the  fluid  reaction,  and  the  cross  spebtral  density  of  surface  pressure  is 

>  !  “  ’’  {  . V  "  i'1' 

5  'V  „  L»  J  :•  X  f! 

0 p  (x,y,to)  =  2ir  <pf  (x,w)p f*  (y,to)> 


The  factor  2  outside  the  summation  accounts  for  the  decorrelation  of  the  upper  and  lower 
surfaces.  The  summation  is  technically  over  all  modes  wljich  are  excited  at  frequency  w. 
Also  in  Equation  (24)  we  have  replaced  Lw  by  As,  denoting  an  area  integration.  The  sur¬ 
face  pressure-cross  spectral  density  depends  on  both  points  of  measurement  x  and  y".  As  is 
discussed  in  the  next  section,  it  may  be  measured  with  pairs  of  pressure  sensors.  In  such  an 

— ■y  ' — ^  — -y  — ^ 

experiment,  the  cross  spectral  density  of  pressures,  measured  at  locations  x  and  y  =  x  +  r , 
is  obtained. 

We  will  simplify  our  discussions  by  expressing  the  pressure  cross  spectrum  in  a  separable 
form  :  i:- 

<Ppp(x;7,w)  =  p2  (x)  $(to)e  e  lrJ/Xz  e  uc  (25) 

In  this  expression  p2(x)  is  the  mean-square  pressure  at  a  chord  wise  coordinate  x,  and  <£(02) 
is  its  autospectrum.  The  pressure  field  is  converted  chordwise  across  the  strut  at  a  speed  U  . 
The  Ax  and  Az  are  each  one-half  of  the  chord-  and  spanwise-integral  correlation  lengths  of 
the  pressure.  The  motivation  for  selecting  this  functional  form  comes  from  the  theories  of 
bound  ary -layer  excitation.  Specifically,  Corcos  (1964)14  has  suggested  the  model  for  the 
pressure  cross  spectrum  as  long  as  the  statistics  of  the  pressure  field  are  homogeneous  in 
space  x.  We  may  expect  Equation  (25)  to  apply  approximately  if  the  field  is  homogeneous 
within  a  correlation  area  2Ax  2Az. 


l^Corcos,  G.M.,  Journal  of  the  Acoustical  Society  of  America,  Vol.  35, ;p.’  192  (1964). 
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In  the  case  of  discrete  vortex  excitation,  an  alternative  form  of  Equation  (25)  is  neces¬ 
sary  to  account  for  the  localization  of  pressure  at  the  trailing  edge.  However,  in  the  experi¬ 
ments  singing  did  not  occur  so  that  case  will  not  be  considered. 

For  the  flow  excitation  of  the  cantilevered  strut,  the  statistics  are  uniform  over  the  span- 
wise  direction  as  long  as  the  turbulent  flow  field  is  unaffected  by  strut  vibration.  The  pres¬ 
sure  field  is  expected  to  change  with  chordwise  distance  x  from  the  leading  edge  because  of 
the  growth  of  the  boundary  layer  on  the  strut.  Near  the  trailing  edge,  where  the  upper  and 
lower  surfaces  converge,  the  static  pressure  gradient  is  adverse,  causing  a  more  pronounced 
thickening  of  the  boundary  layer  in  that  region.  Because  of  the  spatial  nonstationarity  of  the 
flow  field,  we  can  only  approximately  evaluate  Equation  (24)  with  Equation  (25).  For  the 
cases  that  Xx/w  <<  1,  and  X2km  <<  1,  the  integrals  over  rx  and  ry  may  be  approximated, 
giving  for  beam  modes 


0vv(w>zo) 


*m2(z0) 

4\A 

rl/2 

1  Zm  (“)l2 

As 

l/2 

P2  (?) _ .  d  JL 

l+(coXx/Uc)2  w 


(26) 


where  2X  x  and  2Xz  are  the  integral  correlation  lengths  of  the  pressure  field  averaged  over  the 
area  of  the  strut.  If  we  have  measurements  of  the  mean-square  pressure  with  its  correlation 
area  and  autospectrum,  we  can  estimate  the  modal  response  of  the  (m,0)  mode  of  the  beam, 
knowing  the  total  loss  factor  of  the  beam.  The  autospectral  density  for  the  local  acceleration 
of  the  (m,0)  mode  at  resonance  is 


^AA  (“’V  =  2 


<Zo) 


[p  h  +  mm]  277| 


4X^X 

A 


A  r 

A,  J. 


1/2  _ 

P2  (x)  - 

-1/2  l+(“Xx/Uc) 


<J>(co)  X 

- d  — 

,2  W 


(26  B) 


In  Chapter  3  we  will  outline  the  measurement  of  the  statistics  of  the  pressure  field;  in  Chap¬ 
ter  4  we  will  discuss  measurements  of  the  actual  flow-induced— beam-vibration  levels  in  a 

i  , 

water  tunnel  and,  finally,  in  Chapter  5  we  will  compare  those  levels  to  those  which  are  esti¬ 
mated,  using  Equation  (26B)  and  the  statistics  of  the  pressure  field  obtained  in  the  wind 
tunnel. 


3.  WIND  TUNNEL  INVESTIGATION  OF  THE  FLUCTUATING  PRESSURE  FIELD 

A  survey  of  the  turbulent  boundary-layer  mean  and  fluctuating-velocity  profiles  was 
made  in  conjunction  with  two-point,  longitudinal  and  lateral  wall-pressure,  cross  spectral 
densities.  The  measurements  were  performed  in  air  for  zero  and  a  2-deg  angle  of  attack. 
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The  purpose  of  this  investigation  was  to  determine  the  pressure-field  statistics  and  to  ascertain 
which  parameters  may  be  used  in  scaling  the  wind  tunnel  results  to  make  predictions  for  the 
water  tunnel  experiment,  described  in  a  later  section. 


DESCRIPTION  OF  EXPERIMENT 

The  boundary-layer  investigation  was  performed  on  a  strut.  Figure  2a  shows  the  cross- 
section.  The  strut  was  4-ft  in  length  with  a  thickness  of  2  in.  and  chord  of  22  in.  It  was 
positioned  between  two  end  plates  to  insure  two-dimensionality  of  the  mean  flow.  Mean  and 
turbulent  velocity  profiles  were  measured  at  the  lettered  positions  shown  in  Figure  2a.  Auto 
and  cross  spectral  densities  of  wall  pressure  were  measured  at  Positions  A  through  D  on  the 
strut.  A  complete  description  of  the  experimental  technique  and  instrumentation  has  been 
given  by  Blake  and  Dwyer  (1973);15  only  an  outline  will  be  given  here. 

Velocity  measurements  were  made  with  a  Thermosystems  1010A  hot  wire  anemometer, 
used  with  a  Thermosystems  1005B  linearizer.  The  anemometer  probe  was  continuously  tra¬ 
versed  in  directions  normal  to  the  strut.  Incremental  movement  of  the  probe  could  be  made 
in  0.01 -in.  steps.  The  error  in  absolute  positioning  was  generally  about  0.02-in.  because  there 
was  no  satisfactory  means  of  determining  a  positioning  datum  under  flow  conditions  without 
damaging  the  probe.  Using  the  linearized-voltage  outputs,  the  velocity  profiles  could  be  di¬ 
rectly  normalized  on  the  mean  voltage  corresponding  to  the  free-stream  velocity.  The 
boundary-layer  thickness  has  been  defined  as  the  distance  from  the  wall  for  which  the  mean 
linearizer  voltage  is  0.99  of  the  mean  voltage  corresponding  to  the  local  free-stream  velocity. 
Normalized  mean  velocities,  which  are  determined  directly  from  the  linearized  mean  voltages, 
are  in  error  by  less  than  10  percent,  when  U/U*.  «  0.5,  and  by  less  than  20  percent,  when 
U/Uoo  «  0.25.  These  errors  are  basically  determined  by  the  deviation  of  the  relationship  be¬ 
tween  velocity  U  and  voltage  V  from  King’s  law,  which  states 

V2  -  V2  =  Ay/U 


The  anemometer  voltage  at  U  =  0  is  V0,  and  A  is  a  constant.  The  linearizer  adjustments 
were  set  under  the  assumption  that  King’s  law  holds. 

Fluctuating-pressure  measurements  were  made  with  Bruel  and  Kjaer  (B  &  K)  4138  micro¬ 
phones,  which  had  been  fitted  with  pinholes  one-thirty  second  of  an  inch  in  diameter.  The 
pinhole  reduced  the  sensitive  area  of  the  microphone  and  caused  cavity  resonance  at 
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approximately  17  kHz.  The  microphones  were  connected  to  B  &  K  2619  cathode  followers, 
using  the  B  &  K  UA  0123  flexible  right-angle  adaptors,  which  had  been  custom-fitted  with  a 
special  1/8-in.  fitting  for  the  microphone.  The  complete  microphone  system  was  acoustically 
calibrated  in  a  free  field  by  comparison  with  another  microphone  of  known  sensitivity.  The 
calibration  was  reliable  to  ±  1  dB. 

Spectral  analysis  was  performed  online,  using  a  Time  Data  TD  90A  fast  Fourier  proces¬ 
sor.  The  effective  bandwidth  of  analysis  ranged  from  1 2  to  64  Hz  with  total  numbers  of  in¬ 
dependent  samples  ranging  from  256  to  1000.  This  number  of  samples  made  the  statistical 
precision  on  spectral  densities  well  within  ±  10  percent. 


MEAN  AND  TURBULENT 
VELOCITY  PROFILES 

Boundary-layer  surveys  were  made  at  selected  Positions  A  through  G  on  the  strut  and  in 
the  wake  of  the  strut.  Measurements  were  obtained  with  angles  of  attack  a  at  0  and  2  deg. 
When  the  strut  was  set  at  an  angle  of  attack,  the  measurements  were  made  on  the  low  static 
pressure  side  of  the  strut.  Figure  3  shows  normalized  mean-velocity  profiles,  measured  at 
Positions  C  through  F.  The  mean  velocities  are  normalized  on  the  local  free-stream  veloc¬ 
ity;  distance  from  the  wall  is  normalized  on  the  local  boundary -layer  thickness.  The  dotted 
line  in  Figure  3  is  a  mean-velocity  profile,  measured  by  Blake  and  Dwyer15  on  a  similarly 
shaped  strut.  That  strut  was  of  the  same  thickness;  however,  it  had  a  much  longer  parallel 
section  than  the  strut  considered  here.  This  profile  was  measured  in  nearly  zero  pressure  gra¬ 
dient  at  60  ft/sec  and  a  distance  from  the  leading  edge  roughly  corresponding  to  just  forward 
of  Position  C  on  the  strut  considered  here.  Comparison  of  the  zero  pressure-gradient  profile 
with  those  measured  at  Position  C  shows  the  effects  of  the  adverse  pressure  gradient  imposed 
on  the  flow  after  the  knuckle.  Profiles  measured  at  points  downstream  of  Position  C  show 
further  effects  of  the  adverse  pressure  gradient;  the  mean  velocity  near  the  wall  continues  to 
be  retarded;  however,  flow  separation  was  not  observed.  Profiles  measured  at  50  and  100 
ft/sec  inflow  appear  to  be  similar  when  scaled  in  the  manner  of  Figure  3;  however,  profiles 
at  different  streamwise  locations  are  notably  dissimilar  because  of  rapidly  changing  static 
pressure. 

A  flow-visualization  experiment  in  which  the  strut  was  painted  with  a  mixture  of  lin¬ 
seed  oil  and  titanium  dioxide  disclosed  a  small  separated-flow  region  at  the  leading  edge.  The 
region  had  a  streamwise  extent  of  approximately  h/2,  h  being  the  strut  thickness.  The  separa¬ 
tion  bubble  was  centered  on  the  point  at  which  the  circular  leading  edge  was  faired  into  the 
parallel  sides  of  the  strut.  The  flow  visualization  disclosed  no  other  separation  regions  on  the 
strut;  this  was  in  agreement  with  the  anemometer  measurements. 
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Profiles  of  turbulence  intensity,  measured  on  the  strut  at  a  =  0°,  are  shown  in  Figure  4. 

They  generally  reflect  the  effects  of  the  adverse  pressure  gradient  on  the  tapered  section  of 

the  strut.  Figure  4a  shows  a  profile  obtained  in  an  adverse  pressure  gradient  by  Schloemer 

(1967)16  for  qualitative  comparison.  At  all  stations  of  measurement,  as  well  as  in  the 

Schloemer  result,  the  intensities  are  peaked  near  y/5  ~  0.2  to  0.3.  The  measurements  at 

Position  C  were  made  at  Reynolds  numbers  based  on  mean-in-flow  velocity  UQ  and  on  local 

momentum  thickness  N„  =  IL  dlv  on  the  order  of  2.5  x  103;  the  Schloemer  measurement 

was  made  at  Nc  ~  9180.  We  have  used  the  conventional  definition  of  momentum 
Re0 

thickness 


The  large  difference  in  Reynolds  numbers  restricts  the  quantitative  comparison  of  the  Schloemer 
results  with  the  current  ones.  The  results  of  Figure  4  also  show  a  slight,  but  consistent,  re¬ 
duction  in  the  peak-turbulence  intensities  as  the  mean  velocity  is  increased. 

When  the  strut  was  set  at  a  =  2°,  the  boundary-layer  growth  on  the  low-pressure  side  of 
the  strut  was  considerably  altered.  Figure  5  shows  mean-velocity  profiles,  measured  at  Posi¬ 
tions  A,  C,  and  F  at  25  and  50  ft/sec.  Higher  velocities  were  not  used  because  the  strut  was 
a  lifting  surface.  The  slight  angle  of  attack  produced  a  more  severely  adverse  pressure  gra¬ 
dient,  resulting  in  a  much- thickened  boundary  layer.  Furthermore,  the  normalized  velocity 
profiles  are  not  as  well-scaled  using  displacement  thickness  and  free-stream  velocity.  The 
profiles  measured  at  Position  A  are  particularly  dissimilar. 

A  flow-visualization  experiment  at  50  ft/sec  disclosed  an  extensive  separated-flow  region 
at  the  leading  edge  on  the  suction  side  of  the  strut.  The  associated  separation  bubble  appear¬ 
ed  to  extend  about  two  strut  thicknesses  downstream  of  the  leading  edge.  There  were  no 
other  separated-flow  regions  on  the  strut. 

Figure  6  shows  higher  levels  of  turbulence'  intensity  profiles  than  those  at  a  =  0°.  At 
Position  A,  we  note  the  high-intensity  region  extending  to  y/5  ~  0.4;  the  normalized  inten¬ 
sities  are  lower  at  50  ft/sec  than  at  25  ft/sec.  Tins  speed  dependence  is  more  pronounced 
here  than  at  other  positions  at  either  a  =  0°  or  a  =  2°.  The  high  intensities  at  Position  A, 
therefore,  may  be  determined  by  leading-edge  separation,  and  the  pronounced  speed  depend¬ 
ence  of  both  the  mean  and  the  turbulent  intensities  probably  reflects  the  influence  of  speed  on 
the  size  and  strength  of  the  separated  flow.  Further  downstream,  e.g.,  at  Position  C,  normalized 
intensities  are  still  higher  at  a  =  2  °  than  at  a  =  0°;  however,  the  speed  dependence  is  more  similar 
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to  that  observed  at  0  deg.  The  intensities  in  the  near  wake  at  Position  F  were  measured  on  both 
sides  of  the  strut.  The  lack  of  symmetry  reflects  the  effect  of  angle  of  attack.  On  the  pressure 
side  of  the  strut,  the  intensity  profiles  at  0  and  2  deg  are  similar.  The  boundary-layer  thick¬ 
nesses  at  this  location  at  50  ft/sec  are  1.1  in.  for  a  =  0°  and  1.06  in.  on  the  pressure  side  for 
a  =  2°.  On  the  suction  side  at  50  ft/sec  and  a  =  2°,  the  boundary-layer  thickness  is  2.05  in. 
Thus  the  boundary-layer  characteristics  on  the  pressure  side  appear  to  be  quite  unaffected  by 
the  small  change  in  angle  of  attack,  while  the  boundary  layer  on  the  suction  side  is  dramatically 
altered. 

Figure  7  gives  displacement  thickness  5*  measurements.  The  displacement  thickness  is 
calculated  from  the  measured  mean-velocity  profiles  by 


(27) 


We  note  the  generally  higher  values  of.  5*  at  a  =  2°  than  at  a  =  0°.  Table  1  gives  mean- 
velocity  parameters  for  the  strut. 

Wall-shear  coefficients,  defined  as 


Cf  =  TJ* 

where  tw  is  the  mean  local  wall  shear,  were  determined  from  the  mean-velocity  profiles.  The 
mean-velocity  profiles  were  curve  fitted  to  Coles’  (1956)17  law  of  the  wall, 

77-  =  5.75  (28) 

UT 

where  U  is  the  local  mean  velocity, 

UT  is  the  friction  velocity,  and 
v  is  the  kinematic  viscosity  of  air. 

The  friction  velocity  is  related  to  Cf  by  UT/Uoo  =  (Cf/2)^2.  Wall-shear  coefficients  obtained 
in  this  manner  are  subject  to  the  assumption  that  the  law  of  the  wall  applies  in  some  region 
of  the  boundary  layer  close  to  the  wall.  Even  in  severely  adverse  pressure  gradients,  Coles17 
has  shown  Equation  (28)  to  apply.  Some  of  the  mean-velocity  profiles  of  Figures  3  and  5 
have  been  replotted  in  Figure  8.  The  spread  in  the  data  suggests  an  uncertainty  in  calculated 
Cf  values  of  at  least  20  percent.  The  values  are  shown  in  Table  1.  In  the  trailing-edge  region, 
Cf  is  nearly  independent  of  position;  it  is  reduced  by  a  factor  of  1.35  for  a  factor  of  2 
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increase  in  speed.  This  suggests  a  speed  dependence  of  Uq  '4.  At  a  =  2°,  the  speed  depend¬ 
ence  of  Cf  is  less,  and  this  would  imply  a  weaker  Reynolds  number  dependence  than  at 
a  =  0°.  We  note  at  this  point  that  all  turbulence  intensities  peak  at  approximately  y/5  ~ 

0.2  to  0.3.  This  roughly  corresponds  to  y\JT/v  ~  250  for  most  measurement  cases.  The 
mean-velocity  profiles  of  Figure  8  indicate  that  this  coordinate  separates  the  law-of-the-wall 
region  from  the  wake  region  of  the  boundary  layer  where  the  profile  deviates  from  Equation 
(28).  Also,  an  examination  of  Figures  3  and  5  shows  that  the  vertical  gradient  of  mean  veloc¬ 
ity  dtU/U^)  /  d  (y/5)  diminishes  in  most  cases  for  y/5  >  0.3.  Thus  the  high  intensities  are 
restricted  to  the  constant  stress  or  law-of-the-wall  region  of  the  boundary  layer  where  the 
mean  velocity  gradient  is  highest  and  where  the  local  mean  velocity  U  is  less  than  approxi¬ 
mately  0.8  U^. 


STATISTICS  OF  THE  SURFACE 
PRESSURE  FIELD 

Autospectral  densities  of  wall-pressure  fluctuations  were  measured  at  Positions  A  through 
Data  =0°  for  50  and  100  ft/sec.  The  spectra  are  nondimensionalized  on  inflow  speed  UQ, 
dynamic  head  of  the  inflow  q0  =1/2  p0  UQ2  ,  and  5*;  these  variables  collapse  the  spectra 
within  3  dB  as  shown  in  Figure  9.  This  spread  in  the  dimensionless  spectra  is  about  double 
that  shown  by  measurements  in  other  facilities;  however,  the  pressures  measured  on  the  strut 
are  subject  to  a  changing  pressure  gradient.  Thus,  in  the  present  case,  a  unique  dependence 
of  the  spectra  on  dynamic  head  and  displacement  thickness  is  unlikely.  Schloemer16  has 
shown  a  substantial  dependence  of  the  spectra  on  pressure  gradient;  the  effects  of  an  adverse 
pressure  gradient  are  more  severe  than  the  effects  of  a  favorable  one.  In  the  present  case,  we 
have  normalized  the  spectra  at  Position  C  on  the  local  wall  shear  rw  ;  the  right-hand  ordinate 
on  Figure  9  applies  to  this  normalization.  The  dotted  line  in  Figure  9  represents  the  wall- 
pressure  spectra,  normalized  on  r  ,  measured  on  a  wind  tunnel  wall  by  Blake  (1970).18 
Throughout  most  of  the  frequency  range,  the  wind  tunnel  spectrum  is  of  a  level  comparable 
to  the  spectrum  measured  at  Position  C,  and  the  spectrum  shapes  are  nearly  parallel.  At  high 
frequencies,  deviation  of  the  spectrum  shapes  is  caused  by  spatial  averaging  of  the  pressure 
field  over  the  microphone  face.  Corcos14  has  shown  that  if  Uc  is  the  speed  of  convection  of 
the  pressure  field  across  the  microphone,  and  if  Rm  is  the  microphone  radius,  the  pressure 
measurement  is  attenuated  by  roughly  3  dB  at  ojRm  /Uc  ~  1.  We  consider  this  to  be  the 
upper  limit  of  the  nonaveraged  spectrum.  If,  in  anticipation  of  later  results,  we  let  Uc  = 
0.7UQ  ,  and  since  2Rm  «  0.03  in.,  we  find  o>5*/ =  4.6  corresponds  to  the  upper  limit  of 
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the  nonaveraged  spectrum  at  Position  C.  Thus  the  band  limiting  of  the  autospectra  in  Fig¬ 
ure  9  for  a)6*/U0  >  5  is  attributed  to  geometric  microphone-size  averaging  rather  than  a 
nature  of  the  flow. 

The  results  show  that  over  the  downstream  half  of  the  strut,  the  autospectra  are  well 
described  by  the  wall  shear,  8*,  and  UQ,  and  they  are  similar  to  those  measured  in  boundary 
layers  at  higher  Reynolds  numbers  with  and  without  pressure  gradients. 

The  autospectra  of  pressure  measured  at  Position  A  with  a  =  2°  are  shown  in  Figure  10 
for  UQ  =  25  ft/sec,  and  UQ  =  50  ft/sec.  The  autospectra  for  this  position  scale  well  on  speed 
and  displacement  thickness;  however,  the  low-frequency,  dimensionless  spectrum  levels  are 
about  10  dB  higher  than  those  observed  at  all  positions  for  a  =  0°.  Using  the  criterion  that 
was  applied  to  Figure  9,  microphone-size  effects  occur  in  the  spectrum  for  cu8*/U0  >  7. 
Autospectra  measured  at  other  positions  for  a  =  2°  are  shown  in  Figure  1 1 .  The  spectrum 
for  Position  A  is  shown  with  dashes  for  comparison.  The  spread  of  dimensionless  spectra  is 
greater  than  for  a  =  0°,  and  low-frequency  spectrum  levels  appear  to  diminish  as  the  measure¬ 
ment  point  is  moved  further  downstream  and  as  the  flow  speed  is  increased.  It  is  likely  that 
upstream  positions  are  influenced  by  the  leading-edge,  separation  region  which  causes  large- 
scale,  turbulent  eddies  to  be  convected  downstream.  Recall  from  flow  visualization  that  the 
angle  of  attack  increased  the  extent  of  the  leading-edge,  flow  separation.  As  the  flow  be-  \ 
comes  influenced  by  the  converging  surfaces  near  the  trailing  edge,  the  influence  of  the 
leading-edge  separation  is  less  important.  The  spectra  for  co8*/UM  >  1  are  nearly  similar  for 
both  values  of  a,  except  for  possible  microphone  size  attenuation  effects,  which  are  apparent 
at  co5*/U0  >  5,  and  a  =  0°,  and  a>5*/U0  >  7  for  a  =  2. 

Two-point,  cross  spectral  densities  were  measured  with  different  streamwise  separations 
for  both  angles  of  attack.  We  will  adopt  the  representation  of  the  cross  spectrum  originally 
suggested  by  Corcos14  and  will  normalize  the  magnitude  of  the  cross-spectrum  |0(rj,O,co)| 
on  the  auto-spectrum  3>p  (co).  The  magnitudes,  therefore,  represent  a  coherence.  Here  we 
have  let  rj  be  the  streamwise  separation  distance.  The  normalized  magnitudes  will  be  pre¬ 
sented  as  a  function  of  the  measured  phase  7(r2  ,cu),  which  depends  on  both  frequency  and 
separation.  Figure  12  shows  results  at  a  =  0;  the  spectra  were  measured  at  Positions  A  and 
C,  and  the  microphone  separations  were  referred  to  either  of  those  positions.  Upstream 
separations  are  designated  negative;  downstream  separations  are  designated  positive.  Measure¬ 
ments  were  made  at  both  50  and  100  ft/sec.  For  phases  greater  than  three,  the  cross¬ 
spectrum  magnitudes  for  all  speeds  and  separations  are  generally  described  by  the  phase 
alone.  For  phases  less  than  three,  the  magnitudes  are  also  functions  of  the  reference  loca¬ 
tion  of  the  measurement  as  well  as  the  separation.  The  magnitude  of  the  cross  spectrum  is 
slightly  less  than  that  measured  by  Blake18  on  a  smooth  wall  in  zero  pressure  gradient.  The 
lower  magnitudes  are  most  likely  the  result  of  the  adverse  gradient  extant  on  the  strut;  such 
an  effect  has  been  studied  by  Schloemer. 
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The  longitudinal  cross-spectrum  measurements  for  a  =  2°  are  shown  in  Figure  13.  For 
phases  greater  than  two,  the  data  are  scattered;  however,  there  is  no  consistent  dependence  on 
separation  or  measurement  location.  The  magnitudes  are  also  less  than  those  measured  on  the 
smooth,  wind  tunnel  wall  in  zero  pressure  gradient.  At  low  phases,  the  magnitude  of  cross 
spectral  density  is  a  function  of  both  separation  and  the  position  of  reference. 

The  cross-spectrum  phase  has  been  interpreted  by  Corcos14  as 

7j  (rt  ,co)  =  corj  /Uc 

where  Uc  is  the  speed  of  convection  of  a  pressure  component  of  frequency  w  over  the  dis¬ 
tance  rj  .  The  convection  velocities  are  shown  in  Figure  14  for  0-  and  2-deg  angle  of  attack. 
For  o o 5 * / >  0.3,  Uc  /U0  is  nearly  a  constant  independent  of  separation.  The  data  for 
a  2-deg  angle  of  attack  show  slightly  more  scatter  for  t oS*/U„  <  0.75.  At  very  low  fre¬ 
quencies,  Uc  /  U0  >  1  at  Positions  A  and  C.  For  very  low  frequencies  and  large  separations 
at  the  trailing  edge,  Uc  /  UQ  approaches  0.5.  The  convection  velocities  at  high  frequencies 
approach  from  0.6  to  0.7  UQ .  Assuming  that  the  pressure-producing  eddies  are  convected  at 
the  local  mean  velocities;  these  values  of  Uc  place  the  eddies  which  dominate  the  high- 
frequency  spectrum  at  y/5  ~  0.1.  The  eddies  which  dominate  the  low-frequency  spectrum 
are  convected  at  speeds  greater  than  the  inflow  speed  U0 ,  and  they  generate  a  coherent  pres¬ 
sure  field  as  evidenced  by  the  high  values  of  cross  spectral  density  at  low  phase. 

The  measurements  reinforce  the  suggestion  that  the  low-frequency  field  is  generated  by 
the  leading-edge  separation  with  eddies  convected  at  a  velocity  nearly  equal  to  and  sometimes 
larger  than  the  local  free-stream  velocity.  The  low-frequency  pressure  field  is  affected  by 
changes  in  mean  angle  of  attack.  The  high-frequency  pressure  field  is  locally  generated  by 
slowly  convected  eddies  within  the  constant  stress  region  of  the  boundary  layer.  The  pres¬ 
sures  occurring  for  gj5*/U0  >  1  are  only  slightly  affected  by  angle  of  attack. 

The  lateral  cross  spectral  density  measurements  are  normalized  on  the  autospectrum. 
Frequency  is  normalized  on  r3  ,  and  we  take  Uc  ~  0.7  UQ.  The  results  are  shown  in  Figure 
15.  For  phases  ojr3  /UQ  >1,  the  spectral  density  is  similar  to  that  obtained  on  a  wind  tun¬ 
nel  wall  by  Blake.18  At  low  frequencies  the  coherence  is  somewhat  reduced  from  that  ob¬ 
served  in  wind  tunnels,  and  an  additional  dependence  on  separation  appears. 

The  separate  eddy  fields  appear  to  be  characterized  by  separate  coherence  behavior.  At 
higher  frequencies  the  coherence  appears  dependent  on  cor1  /  Uc  and  cor3  /Uc  .  At  lower 
frequencies,  those  apparently  dominated  by  the  separated  flow  at  the  leading  edge,  the  coher¬ 
ence  depends  more  on  the  microphone  separation  than  on  frequency.  Figure  1 6  rather  dra¬ 
matically  illustrates  the  dual  behavior  of  the  cross  spectra.  Frequency  has  been  normalized 
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on  the  strut  thickness  h  and  on  Uo  .  The  ordinate  is  20  log  |  0  (ri ,  0,  w )  I  /  0p  (co)  .  At 
low  frequencies,  the  coherence  for  a  given  separation  is  dramatically  increased  for  angles  of 
attack.  At  a  =  2°,  the  low-frequency  coherence  is  more  a  function  of  separation  than  of 
frequency;  it  is  describable  as 

'  -r./6*/13.7  rl  COh  _  , 

|0  (rx  ,0,co)  1  / <t>p  (co)  ~  e  1  —  <15,  —  <4  (29A) 

0  U0 

Also,  the  lateral  cross  spectral  density  is  described  by 

-r,/5V2.9  r3  coh 

0 (0,r3 ,co) / 0  (co)  ~  e  3  —  <9,  — <4  (29B) 

p  8  uo 

At  high  frequencies,  the  spectra  are  describable  as 

|0  (rj  ,0,co)  l  /4>p  ~e“ai(Wri/Uc)  (30A) 


and 

0(0, r3  ,o j)/0p  (co)  -  e"a3  (wr3  /Uc)  (30B) 

where  a3  and  a3  are  determined  from  curve  fitting  Figures  12,  13,  and  15. 

The  results  of  the  wind  tunnel  investigation  will  be  scaled  to  estimate  the  pressure  field 
developed  on  an  identically  shaped  strut,  studied  in  a  water  tunnel.  Tentatively  based  on  the 
wind  tunnel  results  alone,  the  autospectral  density  of  wall  pressure  appears  to  be  determined 
by  dynamic  head  at  low  frequencies,  which  are  domintaed  by  leading-edge  separation,  and  by 
wall  shear  at  the  higher  frequencies,  which  are  dominated  by  the  local  inner  boundary  layer. 
Displacement  thickness  and  inflow  velocity  appear  to  scale  the  frequency.  The  coherence  of 
the  pressure  field  has  two  separate  descriptions  given  by  Equations  (29)  and  (30). 

4.  WATER  TUNNEL  INVESTIGATION  OF  THE  DAMPING  AND 
VIBRATION  OF  A  BEAM  IN  FLOW 

This  segment  of  the  investigation  was  directed  at  measuring  the  flow-induced  damping 
and  vibration  of  a  series  of  cantilevered  beams  in  water. 
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DESCRIPTION  OF  EXPERIMENT 


Description  of  Test  Specimens 

The  cantilevered  beams  were  generally  20  in.  long  and  were  of  varying  widths.  One  beam 
was  1.5  in.  wide  and  14  in.  long.  The  beams  protruded  from  the  clamp  through  a  6-in.,  dead- 
water  region  into  the  open  jet  of  the  water  tunnel.  A  typical  cross  section  of  the  beams  is 
shown  in  Figure  2b;  the  overall  width  is  designated  w.  The  section  consisted  of  a  circular 
leading  edge  of  0.125-in.  radius,  followed  by  a  1 /4-in. -thick  parallel  section  of  variable  length 
1  ft.  The  trailing  edge  was  a  wedge  with  a  12.5-deg  included  angle.  The  length  of  the  wedged 
section  was  1.0  in.  for  all  beams.  The  widths  of  the  beams  were  varied  by  changing  1  ft;  they 
were  1.5,  2.75,  and  5.0  in.  The  beams  were  stainless  steel. 

Water  Tunnel  and  Damping  Fixture 

The  measurements  were  made  in  the  12-in.,  variable  pressure  water  tunnel  at  the  Naval 
Ship  Research  and  Development  Center.  The  cylindrical  test  section  of  the  tunnel  was  37  in. 
long  and  30  in.  in  diameter  with  a  24-in.-diam  access  hatch,  located  on  the  top  side  of  the 
test  section.  Concentric  with  the  test  section  were  a  9  1 /4-in .-long  entrance  nozzle  and  a 
16-in.-long  exit  nozzle  both  12-in.  in  diameter;  thus  was  produced  a  horizontal  cylindrical  jet, 
approximately  12  in.  in  diameter  and  12  in.  long,  in  which  to  conduct  the  experiment. 

The  clamping  assembly  used  to  hold  the  beams  in  the  waterjet  is  shown  schematically  in 
Figure  17.  Extending  down  13  3/4  in.  from  the  top  access  hatch  cover  were  four  1  1/4-in. 
angle-iron  brackets  at  the  bottom  of  which  were  attached  two  1  1/4-  by  13  3/4-in.  angle-iron, 
horizontal  cross  members.  A  190-lb,  mass-clamp  system —consisting  of  a  10-  by  10-  by  7-in. 
stainless  steel  block  weighing  185  lb  and  a  2-  by  5-  by  5-in.  clamping  block  weighing  5  lb, 
similar  in  design  to  that  used  by  Kerlin  and  Snowdon  ( 1 97 0)1 9  —  was  mounted  on  the  two 
cross  members  by  means  of  three  cylindrical  vibration  mounts  from  Barry  Controls,  Type 
A32-051.  Tire  system  had  resonant  frequencies  of  12  Hz  ±42  percent  in  the  vertical  direction 
and  5  Hz  ±  16  percent  in  the  horizontal  direction.  These  frequencies  were  estimated,  using 
stiffnesses  provided  in  the  data  sheet  supplied  with  the  vibration  mounts. 

Extending  down  vertically  from  the  clamping  block,  the  beams  were  clamped  rigidly  with 
a  root  depth  of  5  in.  between  the  blocks,  using  six  3/8-16  cap  screws  tightened  to  40  ft  lb. 

The  cap  screws  were  always  tightened  to  the  same  torque  in  order  to  insure,  as  nearly  as  pos¬ 
sible,  the  same  boundary  condition  for  all  beam  specimens.  Mounted  6  in.  from  the  clamped 
end  of  the  beam  was  a  0.38-in.-diam,  0.5-in.  long,  7-g  Wilcoxon  Research  Model  165  acceler¬ 
ometer  with  a  sensitivity  of  9  mV/g  ±  1  dB. 


l^Kcrlin,  R.L.  and  J.C.,  Snowdon,  Journal  of  the  Acoustical  Society  of  America,  Vol.  47,  p.  222  (1970). 
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Vibration  Analysis 

Accelerometer  signals  were  carried  by  low-noise  coaxial  cable,  which  passed  through  a 
packing-gland  seal  in  the  test-section  wall.  The  signals  were  then  amplified  by  an  Ithaco 
Model  257-A  preamplifier  and  were  processed  in  either  of  two  ways.  On  one  hand,  spectral 
densities  of  acceleration  were  obtained,  using  a  General  Radio  (GR)  Model  1900-A  wave  ana¬ 
lyzer  with  a  3-Hz-bandwidth  filter.  Spectra  were  recorded  on  a  GR  Model  1521-B  graphic 
level  recorder.  On  the  other  hand,  damping  measurements  were  obtained  by  filtering  decay¬ 
ing  acceleration  levels,  using  a  Bruel  and  Kjaer  (B  &  K)  audiofrequency  spectrometer  Type 
2111,  and  by  recording  the  filtered  levels  on  a  B  &  K  Type  2305  graphic  level  recorder.  Loss 
factors  for  modes  were  determined  from  the  1 /3-octave-band,  60-dB  decay  times  resulting  from 
the  response  of  the  beam  to  the  impulses  received  from  a  mechanical  impulse  generator.  The 
structural  mode  densities  of  the  beams  were  small  enough  that  a  1 /3-octave-measurement 
band,  centered  on  a  given  mode,  did  not  include  more  than  a  single  beam  mode. 

Mode  shapes  of  the  beams  in  still  water  were  measured,  using  a  noncontact  displacement 
probe  called  a  Fontonic  Sensor,  Model  KD45A,  Mechanical  Technology,  Inc. 


Anemometer  Positioning  and  Flow  Measurement 

One  of  the  two  17-in.-diam  viewing  portholes  on  the  side  of  the  test  section  was  modi¬ 
fied  so  that  a  7/8-in.  brass  tube,  which  held  a  DISA  55A83  hot  film  anemometer  probe,  could 
be  inserted  into  the  tunnel  normal  to  the  beam  surface.  The  viewing  port  was  1  l/4-in.-thick 
Plexiglas  and  had  a  4-in.-diam  penetration  hole,  which  was  covered  with  a  6  3/4-in.-diam,  1/8- 
in.-thick,  brass  plate  that  held  the  7/8-in.  tube.  The  tube  was  held  in  a  combination  knurled 
screwdown  clamp  and  0-ring  seal  which  was  brazed  to  the  plate.  The  brass  plate  was  sealed 
against  the  Plexiglas  via  an  0-ring  and  was  held  against  the  plexiglas  with  three  clamp  screws. 
Likewise,  the  anemometer  probe  and  the  7/8-in.  brass  tube  were  sealed  with  0-rings.  The 
brass  plate-hole  cover  enabled  the  anemometer  to  be  moved  along  the  directions  of  the  beam 
chord  and  length,  while  the  0-ring  seal  around  the  7/8-in.  brass  tube  allowed  the  anemometer 
to  be  positioned  along  the  normal  to  the  beam.  The  position  measurements  were  made  with 
a  machinist  rule  with  0.01 -in.  graduations.  A  DISA  55D05  anemometer  was  used  to  drive  the 
anemometer  probe. 

The  mean  inflow  velocity  was  measured  with  a  pitot-static  tube  and  a  mercury-in-glass 
manometer.  This  velocity  measurement  provided  the  anemometer  calibration  standard.  The 
probe  was  mounted  to  the  frame  of  the  clamping  block  so  that  the  probe  tip  was  just  up¬ 
stream  of  the  leading  edge  of  the  beam. 
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Mechanical  Impulse  Generator 


In  air  a  mechanical  impulse  could  be  delivered  by  striking  the  beam  with  a  hard,  hand¬ 
held  object.  This  method,  of  course,  was  impossible  in  the  closed  and  water-filled  test  cham¬ 
ber.  The  first  design  of  a  remotely  powered,  impulse  generator  consisted  of  a  solenoid  sur¬ 
rounding  a  cylindrical  plunger.  The  solenoid  was  in  a  series  circuit  with  a  battery  and  micro¬ 
switch.  This  arrangement  produced  impulses  which  were  too  long  because  of  the  long  time 
duration  of  the  pulse  current.  The  plunger  of  the  solenoid  in  this  arrangement  would  chatter 
against  the  beam,  regardless  of  the  speed  with  which  the  microswitch  was  activated.  Hence, 
it  was  decided  to  discharge  a  capacitor  through  a  solenoid  in  order  to  produce  a  short  time 
current  pulse.  A  small  RC  time  constant  was  needed  to  produce  a  pulse  of  short  time  dur¬ 
ation,  while  a  large  capacitance  was  needed  to  provide  enough  energy  at  a  given  charge  volt¬ 
age  to  excite  the  plunger  during  discharge.  Consequently,  it  was  necessary  to  find  a  solenoid 
with  a  resistance  as  low  as  possible,  capable  of  carrying  a  large  current  lead.  Fortunately,  a 
starting  solenoid  from  an  old  Chevrolet  automobile  was  available,  and  this  was  modified  to 
produce  the  mechanical  impulse  generator.  The  two  separate  windings  of  the  solenoid  were 
reconnected  in  series;  then  the  entire  solenoid  along  with  the  lead-in  cable  was  waterproofed. 

The  plunger  consisted  of  a  brass  rod  1/4  in.  in  diameter  and  10  in.  in  length,  which  was 
held  concentric  with  an  iron  pipe  3/4-in.  in  diameter  and  1/4  in.  in  length  by  using  a  brass 
bushing.  The  brass  rod  protruded  from  the  end  of  the  solenoid.  Figure  18  shows  a  sketch 
of  the  assembly.  Parallel  to  the  axis  of  the  rod  the  cross  bushing  was  drilled  out  as  much  as  pos¬ 
sible  to  reduce  both  the  mass  of  the  plunger  and  the  frictional  force  of  the  water  on  it.  The 
plunger  was  pulled  into  the  solenoid  by  the  magnetic  field,  and  a  spring  was  used  to  keep  the 
plunger  away  from  the  stop  of  the  solenoid  after  the  magnetic  field  had  collapsed.  Power  for  the 
solenoid  was  provided  by  a 4000  /rF,  50-V  capacitor,  which  was  charged  to  the  desired  poten¬ 
tial,  using  a  Sorensen  Model  QRD40-0.75  power  supply.  Finally,  a  double  pole,  double  throw 
switch  was  used  to  isolate  the  power  supply  from  the  discharge  transient.  The  final  assembly 
had  a  resistance  of  only  1.5  ohm. 

In  operation,  the  capacitor  was  charged  to  the  desired  potential  then  discharged  via  the 
double  pole,  double  throw  switch,  through  the  solenoid,  producing  a  strong,  short  time  dura¬ 
tion  magnetic  field,  which  pulled  the  plunger  into  the  solenoid  and,  consequently,  caused  the 
tip  of  the  1/4-in.  brass  rod  to  strike  the  beam.  The  strength  of  impact  was  changed  by  vary¬ 
ing  the  charge  on  the  capacitor. 

Test  Procedure 

Before  a  typical  run,  the  water  tunnel  was  filled  and  then  deaerated  for  1  hr  to  remove 
bubbles  from  the  test  chamber.  Then,  following  selection  of  a  jet  velocity,  the  3-Hz-bandwidth 
levels  of  flow-induced  beam  acceleration  were  obtained.  Finally  the  impulse-excited  decay  times 
were  6  measured. 


RESULTS  OF  MEASUREMENTS 


Resonance  Frequencies  and  Mode  Shapes 

The  fundamental  beam  mode  m  =  0  for  a  20-in.-long  by  0.2-in. -thick  cantilevered  beam 
has  a  theoretical  resonance  frequency  in  air  of  approximately  14  Hz;  however,  the  actual  fun¬ 
damental  frequency  could  not  be  measured  by  shaking  because  that  mode  was  outside  the 
frequency  limitations  of  the  instrumentation.  The  resonance  frequency  for  the  2.75-in.  beam 
in  water,  calculated  using  Equations  (9)  for  the  entrained  mass,  is  about  10  Hz.  Narrow 
band,  filtered  levels  of  flow-induced  acceleration  on  some  beams  displayed  peaks  in  the  range 
from  10  to  15  Hz  which  became  less  distinct  as  water  speed  was  increased.  Undoubtedly,  the 
fundamental  resonance  was  influenced  by  the  finite  mass  of  the  clamp  as  well  as  by  the  stiff¬ 
ness  of  the  rubber  mounts.  The  theoretical  resonance  frequencies  in  water  are  given  by  a 
combination  of  the  relationships  following  Equations  (5)  and  (17),  i.e. 

-  (2m  +  l)2  ,  m  >  1  (31) 

4L2 

Measured  resonance  frequencies  are  given  in  Table  2;  for  beam  (m,0)  modes  in  air,  they  are 
higher  than  those  calculated  by  using  Equation  (31)  by  less  than  15  percent.  Torsion  mode 
(0,n)  resonance  frequencies  were  not  calculated.  Entrained  mass  for  the  bending  modes,  de¬ 
termined  from  the  resonance  frequencies  of  Table  2,  was  in  agreement  with  Equation  (9). 

Nodal  patterns  in  water  for  beam  modes  of  the  2.75-  by  20-in.  beam  are  shown  in  Fig¬ 
ure  19.  Locations  of  the  shaker  and  reference  accelerometer  are  designated  by  S  and  A, 
respectively,  on  each  diagram.  Also,  the  node  lines,  which  are  the  theoretical  points  calcu¬ 
lated  by  Young  and  Felgar  (1949), 20  are  shown  for  each  resonance.  Generally,  the  measured 
nodal  patterns  are  curved,  the  curvature  is  caused  by  both  the  varying  thickness  and  the  fi¬ 
nite  width  of  the  beam.  For  the  m  =  2  mode,  the  mode  shape  in  water  was  measured  for 
two  chordwise  locations  along  the  knuckle  and  near  the  center  of  the  beam.  The  arrows  in 
Figure  19  show  the  measurement  locations.  The  results  are  compared  with  the  theoretical  in 
vacuo  mode  shape  of  Young  and  Felgar  in  Figure  20.  The  vibration  levels  were  normalized 
to  be  two  at  the  tip.  With  this  normalization,  the  curves  of  Figure  20  are  consistent  with 
the  normalized  functions  of  Equation  (5).  We  see  that  the  theoretical  mode  shape  for  m  =  2 
is  only  a  rough  approximation  to  the  actual  vibration  patterns  at  the  resonance  frequency. 
The  lack  of  similarity  is  probably  due  to  the  variable  thickness  of  the  beam  and  to  the 

20young,  D.  and  R.P.  Felgar,  “Tables  of  Characteristic  Functions  Representing  Normal  Modes  of  Vibration  of  a  Beam,” 
University  of  Texas,  Publication  4913,  Engineering  Research  Series  44  (Jul  1,  1949). 
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non-ideal  clamping,  since  measured  vibration  patterns  in  air  are  similar  to  those  measured  in 
water.  Theoretical  in-vacuo  mode  shapes  for  the  other  modes,  obtained  from  Reference  20, 
are  also  shown  in  Figure  20.  The  single  point  acceleration  measurement,  located  on  all  beams 
at  z/L  ~  0.28  and  denoted  by  arrows  in  Figure  20,  can  be  considered  representative  of  the 
theoretical  mean-square  level.  For  the  m  =  4  mode,  the  deviation  of  the  mode  shape  from 
the  theoretical  nodal  pattern  as  shown  in  Figure  19  probably  accounts  for  the  acceleration 
levels  being  measured  at  this  location  for  this  mode.  Nevertheless,  measured  vibration  levels 
for  the  m  =  4  mode  are  likely  to  be  relatively  low.  Also,  since  the  beam  passed  through  a 
dead-water  region  into  the  jet,  the  flow  covered  z/L  >  0.3  for  the  20-in.-long  beams  and  z/L 
>  0.43  for  the  14-in.-long  beam. 

Figure  21  shows  3-Hz-bandwidth  acceleration  levels  for  the  2.75-in.  beam  at  23.5  ft/sec 
with  the  mode  orders  responsible  for  the  peaks  shown  in  parentheses.  The  frequency  is  on  a 
linear  scale,  and  the  acceleration  levels  are  referred  to  the  gravitational  acceleration.  The  ac¬ 
celeration  peaks  at  low  frequencies  have  larger  band  widths  than  at  high  frequencies  because 
of  the  higher  loss  factors  at  low  frequencies.  Other  spectra  at  low  speeds  show  a  peak  at 
approximately  10  to  15  Hz,  which  corresponds  to  the  fundamental  mode  of  motion.  At  this 
point  it  is  well  to  note  the  coincidence  of  the  expected  fundamental  with  approximate  reso¬ 
nance  frequency  of  the  rocking  mode  of  the  block  1 2  Hz  and  the  1 2-Hz  rotation  rate  of  the 
tunnel  drive-propeller  blades  at  UQ  ~  7  ft/sec.  Thus  we  had  at  nearly  7  ft/sec  the  well- 
defined  1 0-  to  1 5-Hz  peak  in  the  acceleration  spectrum  which  could  account  for  a  ±  0.5-deg 
angle  of  attack  oscillation. 


Dependence  of  Damping  on  Flow 

Measured  total  loss  factors  for  all  of  the  20-in.  beams  in  still  water  are  shown  in  Figure  22. 
The  loss  factors  show  no  consistent  dependence  on  the  width  of  the  beam,  and  they  are  high¬ 
er  in  water  than  in  air,  where  they  are  less  than  1 0“ 3 .  The  effect  on  damping  of  wetting  the 
clamp  was  determined  in  a  small  tank  with  and  without  immersion  of  the  clamp.  Only 
enough  water  was  used  to  cover  the  clamp,  leaving  the  complete  length  of  the  beam  exposed 
to  air.  The  loss  factors  were  unaffected  by  wetting  the  clamp.  Thus,  it  is  assumed  that  the 
damping  of  the  beam  in  air  is  controlled  by  dry  friction  in  the  clamp.  The  damping  caused 
by  viscous  losses  is  estimated,  using  an  empirical  result  obtained  on  blunt-edged,  free-free 
beams  by  Blake.1  Designating  the  viscous  loss  factor  by  tjv,  we  have  from  Reference  1 
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This  expression  has  been  plotted  in  Figure  22,  using  parameters  for  the  2.75-in.-wide  beam. 
The  observed  loss  factors  are  at  least  three  times  the  calculated  ones.  It  is  possible  that  the 
sharp  trailing  edge  has  increased  the  viscous  losses.  Most  of  the  results  shown  in  Reference  1 
were  obtained  with  blunt-edged  beams,  and  the  effect  of  edge  sharpness  was  never  fully  ex¬ 
amined.  In  an  experiment  comparing  the  round-  to  sharp-edged  damping  of  a  cantilevered 
beam,  Kerr,  Shannon,  and  Arnold  (1940)21  determined  a  factor  of  3.5  increase  in  damping 
level  with  sharp  edges.  In  the  current  experiments,  higher  damping  was  observed  when  small 
bubbles  were  attached  to  the  beam  surface.  These  were  generally  removed  after  deaeration 
of  the  tunnel.  Figure  23  shows  the  increase  of  loss  factor  with  inflow  speed  U0  for  the  2.75- 
in.  beam.  The  damping  is  highest  at  low  frequencies,  approaching  nearly  9  percent  of  critical 
for  the  m  =  1  mode. 

The  data  have  been  normalized  in  line  with  the  interpretation  of  Equations  (18)  and 
(19).  The  measured  loss  factors  are  multiplied  by  the  ratio  of  the  fluid  plus  structural  mass¬ 
es  to  the  entrained  fluid  mass.  At  high  enough  speeds  when  flow  effects  dominate  the  damp¬ 
ing,  we  can  interpret  measurements  as  hydrodynamic  loss  factors 


which  are  given  by  Equation  (19).  Frequency  is  nondimensionalized,  using  the  inflow  veloc¬ 
ity  UQ  and  the  beam  width  w.  Hydrodynamic  loss  factors  for  both  the  1.5-  and  2.75-in.- 
wide  by  20-in.-long  beams  are  shown  in  Figure  24.  The  data  from  Figure  21  are  included 
and  cover  bending  modes  m  =  1  through  m  =  4  as  well  as  the  first  torsion  mode  of  the  2.75- 
in.  beam.  At  very  low  reduced  frequencies-high  speeds  and  low  frequencies-the  loss  factors 
were  calculated  from  the  bandwidths  of  the  flow-induced  acceleration  spectra.  It  appears 
that  the  hydrodynamic  loss  factor  is  well  described  as  a  function  of  reduced  frequency 
ww/U0  alone.  Also  shown  in  Figure  24  are  the  two-  and  three-dimensional  asymptotic  ap¬ 
proximations  for  the  hydrodynamic  loss  factors  given  by  Equations  (20)  and  (21).  Since 
only  a  fraction  of  the  length  of  the  beam  1  /L  is  exposed  to  the  waterflow,  calculated  r?H 
includes  this  factor.  At  very  low  reduced  frequencies,  the  data  approach  the  two-dimensional 
theoretical  loss  factors,  and  the  measured  loss  factors  approach  the  theory  for  vanishing  as¬ 
pect  ratio  /R  at  high  reduced  frequencies.  At  still  higher  reduced  frequencies,  those  consist¬ 
ent  with  small  values  of  UQ ,  the  measured  values  of  rjH  approach  an  asymptotically  con¬ 
stant  value  equal  to  the  Stillwater  loss  factors. 


21}Cerr,  W.  et  al.,  “The  Problems  of  the  Singing  Propeller,”  Proceedings  of  the  Institute  of  Mechanical  Engineering  Vol. 
144,  p.  54  (1940). 
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Figure  25  shows  hydrodynamic  loss  factors  for  the  5-in.  wide  by  20-in.-long  beam.  In 
this  case  data  for  beam  modes  m  =  1,  m  =  2,  and  m  =  4  as  well  as  for  the  torsion  mode  (0,1) 
and  plate  modes  (1,1)  and  (2,1)  are  shown.  Theoretical  values  given  by  Equations  (20)  and 
(21)  are  included  for  the  beam  modes;  the  theoretical  two-dimensional  loss  factor  for  the 
torsion  (0,1)  mode  given  by  Equation  (22)  is  also  included.  The  loss  factors  for  both  mode 
types  approach  the  two-dimensional  values  at  low  reduced  frequencies.  Figure  25  shows  the 
hydrodynamic  loss  factors  for  the  1.5-  by  14-in.  beam.  Only  beam  modes  were  examined, 
and  the  loss  factors  are  again  dependent  only  on  reduced  frequency.  Also  included  are  the 
theoretical  two-  and  three-dimensional  asymptotic  limits  as  well  as  a  more  precise  computa¬ 
tion  using  results  tabulated  by  Lawrence  and  Gerber.11  The  excellent  agreement  with  theory 
as  well  as  the  close  collapse  of  the  data  for  all  beams  is  strong  evidence  for  the  simple  de¬ 
scription  of  the  hydrodynamic  loss  factor  as  a  function  of  reduced  frequency  alone. 

The  loss  factors  at  high  reduced  frequencies  in  Figures  24  through  26  are  influenced  by 
the  damping  extant  in  the  beam  at  zero  flow.  We  could  delete  the  contribution  at  zero  flow 
by  assuming  that  the  total  loss  factor  is  a  linear  sum  of  that  hydrodynamically  induced  and 
that  extant  at  UQ  =  0.  Unfortunately  such  a  correction,  although  appealing,  is  poorly  found¬ 
ed  since  we  have  not  characterized  the  mechanisms  of  damping  in  water  at  U0  =  0.  Thus  we 
cannot  say  that  the  same  mechanisms  exist  unchanged,  when  UQ  f  0.  Viscous  damping,  for 
example,  may  well  be  modified  when  a  steady  flow  is  imposed  on  the  beam. 

Another  effect,  until  now  overlooked,  is  that  caused  by  finite  thickness.  The  two-  and 
three-dimensional  unsteady  foil  theory  has  been  developed  for  struts  of  vanishing  thickness. 
When,  in  fact,  the  reduced  frequency  is  high  enough  so  that  cuh/UQ  >  1 ,  we  can  expect 
thickness  effects  to,  perhaps,  change  the  observed  damping.  These  points  correspond  to 
(  ww/U0  >  20,  11,  and  6  for  the  5-,  2.75-  and  1.5-in.  beams,  respectively.  If  we  subtract 
loss  factors  measured  at  U0  =0  from  those  measured  at  UQ  f  0  we  obtain  values  that  are 
considerably  lower  than  the  theoretical  lines  shown  in  the  figures  for  cow/U0  >  40  or  so. 

This  has  two  possible  causes:  thickness  effects  reduce  the  hydrodynamic  damping,  and  the 
loss  factors  at  U0  =  0  do  not  apply  at  UQ  f  0  because  viscous  effects  on  the  damping  are 
reduced,  thus  lowering  the  residual  damping. 


Flow-Induced  Vibration  Measurements 

The  records  of  3-Hz-band  acceleration  levels,  similar  to  the  levels  shown  in  Figure  20, 
were  analyzed  to  determine  modal  acceleration  spectral  densities  at  each  resonance  frequen¬ 
cy.  These  spectral  densities  $AA  (co,  xQ)  were  then  nondimensionalized  on  flow  variables  to 

m 

emphasize  the  functionality  demonstrated  by  Equation  (26B).  Letting  the  mean-square,  3-Hz- 
band  levels  measured  at  a  position  on  the  beam  xQ  be  a2,  the  spectral  density  is 


.28 


where 


or 


Af  =  i?xf„ 


Af  =  3 


t?t  f  <3  Hz 

1 1  m 


1?T  fm  >3  HZ 


where  rjt  is  the  total  loss  factor  of  the  mode  m,  and  fm  is  its  resonance  frequency.  The  re¬ 
sulting  acceleration  spectra  a2  /Af  are  shown  for  the  2.75-in.  beam  in  Figure  27.  With  the 
exception  of  some  anomalous  behavior  at  the  lowest  speed,  the  spectra  vary  smoothly  with 
frequency.  We  note  that  only  beam  (m,0)  mode  response  was  examined  in  this  manner. 
Mode  orders  are  shown  in  parentheses  in  Figure  27. 

A  rearrangement  of  Equation  (26B)  to  emphasize  nondimensionalization  is 
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where  we  have  replaced  pph  by  the  area  density  mb  and  have  introduced  the  dynamic  head 
of  the  inflow.  Frequency  should  be  nondimensionalized  on  a  boundary-layer  length  scale  as 
emphasized  in  Chapter  3;  however,  the  magnitude  of  that  scale  is  unknown.  Figure  7  indi¬ 
cates  that  the  dominant  boundary-layer  growth  occurs  at  the  wedge  of  the  cross  section  so 
that  we  have  used  the  maximum  thickness  of  the  section  h  as  a  characteristic  length  for 
scaling  purposes.  The  dimensionless  spectra  for  the  2.75-in.  beam,  taken  from  Figure  27,  are 
shown  in  Figure  28.  The  flagged  points  denote  data  for  the  n  =  4  mode.  The  nondimension¬ 
alization  of  Equation  (32)  collapses  the  data  for  individual  speeds  very  well.  Close  corre¬ 
spondence  of  the  points  suggests  that  the  boundary-layer  length  scale,  which  we  have  previ¬ 
ously  assumed  proportional  to  h,  and  Xx  and  \  are  only  weak  functions  of  speed  and  are 
not  dominantly  affected  by  the  vibration  level.  The  data  for  the  n  =  4  mode  may  be  rela¬ 
tively  low  because  of  the  proximity  of  the  accelerometer  to  the  theoretical  and  actual  loca¬ 
tions  of  the  node  line  of  that  mode. 


The  nondimensionalization  was  repeated  for  the  other  beams  of  the  set;  Figure  29  shows 
results  for  the  5-in.  beam.  Although  the  points  for  the  5-in.  beam  collapse  well  in  terms  of 
U0  and  h,  the  dimensionless  spectra  deviate  slightly  from  those  of  the  2.75-in.  beam  at  high 
frequencies.  This  would  indicate  that  the  actual  length  scale,  which  determines  the  high- 
frequency  spectrum,  is  about  20  percent  smaller  relative  to  h  on  the  5-in.  beam  than  it  is  on 
the  2.75-in.  beam.  Figure  30  shows  dimensionless  spectra  for  the  1.5-  by  20-in.  beam;  the 
spread  in  the  data  includes  the  mean  lines  for  both  the  other  beams.  In  the  case  of  the  1 .5- 
by  20-in.  beam,  the  m  =  4  node  line  is  apparently  quite  near  the  accelerometer  since  the 
flagged  points  are  consistently  the  lowest.  Finally,  Figure  31  shows  dimensionless  spectra  for 
the  1.5-  by  14-in.  beam.  The  dimensionless  accelerations  are  all  in  agreement.  The  mean  line 
for  the  2.75-in.  beam  has  been  linearly  adjusted  for  the  difference  in  flow  coverage  over  the 
span  denoted  by  lw  . 

All  the  normalized  modal  acceleration  spectral  densities  appear  to  be  nearly  unique  func¬ 
tions  of  coh/UQ  alone.  To  examine  the  possible  analytical  modeling  of  the  acceleration  spec¬ 
trum,  we  require  a  combination  of  Equations  (29)  and  (30)  with  Equation  (32),  using  the 
appropriate  fluid-scaling  laws. 


Mean  and  Turbulent  Velocity  Profiles  Behind  the  Beam 

To  establish  a  basis  of  scaling  the  flow  characteristics  obtained  in  air  to  the  conditions 
existing  in  the  water  tunnel,  it  was  necessary  to  make  a  cursory  velocity  survey  in  the  vicinity 
of  the  flow-excited  beam.  The  measurement  was  performed  for  two  speeds  at  Position  F  of 
the  2.75-in.  beam  at  a  distance  z/L  —  0.8  below  the  clamping  block. 

The  mean-velocity  profiles  were  measured  at  two  flow  speeds  of  3.7  and  7.7  ft/sec.  The 
data  were  obtained  on  both  sides  of  the  beam.  Figure  32  shows  all  four  profiles  normalized 
as  U/U^  versus  y/5,  where  y  is  measured  from  the  center  of  the  axis.  The  profiles  agree 
with  each  other  and  with  the  profiles  measured  in  air.  The  strut  was  about  1-deg  constant 
angle  of  attack,  and  the  very  low  frequency  vibratory  motion  of  the  beam  accounted  for 
another  ±0.5  deg  superimposed  on  this.  Thus  the  instantaneous  angle  of  attack  was  bounded 
between  a  0-  and  2-deg  angle  of  attack. 

Figure  33  shows  longitudinal  turbulence  intensities  at  Position  F  behind  the  2.75-in. 
beam.  Included  are  the  results  from  Figures  4  and  6  for  Position  F  obtained  in  air;  the  re¬ 
sults  obtained  in  water  are  generally  bounded  by  the  0-  and  2-deg  angle-of-attack  results  in 
air.  At  values  of  y/5  >2  the  turbulence  intensity  is  approximately  0.01.  This  value  was 
obtained  with  a  5-Hz,  high-pass  filter;  with  a  100-Hz,  highpass  filter,  the  background  inten¬ 
sity  is  only  0.002.  The  autospectral  densities  of  the  turbulence  levels  obtained  at  Position 
F  in  air  and  in  water  are  shown  in  Figure  34.  The  measurements  were  obtained  well  within 
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the  boundary  layers.  The  autospectra  are  normalized  on  the  local  mean  velocity  U,  local 
mean-square  velocity  u2  ,  and  8*.  The  normalized  spectra  display  an  overall  spread  of  6  dB. 
All  measurements  at  y  =  0  are  in  agreement,  regardless  of  the  flow  medium  or  angle  of  attack. 
At  y  —  28*  the  measurements  in  air  are  slightly  higher  at  midfrequencies.  However,  they 
coalesce  better  at  low  and  high  frequencies.  The  close  similarity  of  turbulent  velocities  in  air 
and  water  further  substantiates  the  flow-similarity  in  both  media. 

Table  3  summarizes  some  of  the  boundary-layer  characteristics  at  Position  F  in  air  for 
«  =  2°  and  in  water;  8/8*  and  8*/h  are  nearly  invariant  over  the  Reynolds  number  range  in¬ 
dicated.  The  experiment  was  designed  to  maintain  similar  air  and  water  Reynolds  numbers 
so  that  the  scaling  of  results  would  be  simple. 

The  background  turbulence  spectra  were  obtained  for  frequencies  above  5  Hz  at  three 
speeds.  In  Figure  35  the  spectral  densities  are  shown  nondimensionalized  as 
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where  h  is  a  length  arbitrarily  selected  to  be  equal  to  the  beam  thickness  and  h  =  0.25  in. 

For  speeds  of  UQ  =  7  ft/sec  and  U0  =  13.2  ft/sec  this  nondimensionalization  seems  to  gener¬ 
ally  describe  the  spectra;  for  UQ  =  3.5  ft/sec  this  spectrum  level  is  uniformly  8  dB  lower  than 
those  at  the  higher  speeds.  For  later  use  we  will  determine  a  longitudinal  integral  scale  by  assum¬ 
ing  that  Taylor’s  Hypothesis  of  frozen  convection  applies  to  the  turbulence  and  that  the  spectrum 
is  of  the  form 
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The  background  intensity  is  u£  and  L  is  an  integral  scale  of  the  intensity.  Matching  this 


spectrum  form  to  the  high-speed  data  in  Figure  35  and  using  u| 
L'  -  2h. 


-tW2 


/Un  =  0.03  we  find 


These  spectra  are  compared  to  those  determined  in  the  wake  of  the  2.75  in.  beam  as 

shown  in  Figure  34.  We  multiply  the  high-speed  data  of  Figure  35  by  u2/ug  where  u2  is 

here  taken  as  the  typical  mean-square  intensity  in  the  wake  and  which  is  on  the  order  of 

—1/2 

0.01  U2  .  We  will  also  take  Ug  /UQ  =  0.028;  this  is  the  broadband  intensity  at  tunnel  speeds 
of  7  ft/sec  and  13  ft/sec.  We  will  replace  h  by  8*  =  0.05  in.  So  modified,  the  high-speed 
spectrum  of  Figure  35  has  been  plotted  in  Figure  34.  This  background  is  not  influenced  by 
electronic  noise  and  we  see  that  except  at  the  lowest  frequencies  the  shear  layer  turbulence 
is  not  comparable  to  the  water  tunnel  background. 
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5.  ANALYTICAL  ESTIMATES  OF  THE  FLOW-INDUCED  ACCELERATION  SPECTRUM 

RESPONSE  TO  THE  BOUNDARY-LAYER  PRESSURE 

An  approximate  expression  for  the  acceleration  induced  by  the  boundary  layer  generated 
on  the  strut  can  be  obtained  from  Equation  (32).  We  let  i//m2  (zQ)  1,  and  rewrite  Equation 
(32),  using  the  nomenclature  of  Figures  9  through  1 1 
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where  ?.  =  1  or  2,  depending  on  whether  the  flow  is  correlated  on  one  side  or  on  both  sides 
of  the  beam.  In  the  case  of  the  nonzero  angle  of  attack,  the  low-frequency,  boundary-layer 
pressures  are  mostly  generated  on  one  side,  while  those  at  high  frequencies  are  equal  and  un¬ 
correlated  on  both  sides.  Tire  area  of  the  beam  which  supports  the  most  intense  pressure 
field  is  denoted  by  ww  £w  ,  where  £w  is  the  depth  of  immergence  into  the  wateijet  of  the 
tunnel.  The  close  similarity  of  the  flow  fields  at  Position  F  in  air  and  in  water  as  indicated 
in  Table  3  make  scaling  the  results  of  Chapter  2  straightforward.  Since  the  beam  is  at  an 
angle  of  attack,  the  data  of  Chapter  2  for  a  =  2°  will  be  used.  The  observation  that  at  Posi¬ 
tion  F  5*/h  varies  only  slightly  over  the  Reynolds  number  range  of  Table  3  indicates  that  the 
average  of  8*  over  the  wedge  portion  of  the  beam,  say  5*  should  be  nearly  proportional  to 
h.  At  25  ft/sec,  NRe  t  =  2.9  x  105  in  the  air  experiment,  and  NRe  =  1.6  x  105  in  the  water 
experiment  at  7.7  ft/sec.  Over  this  Reynolds  number  range,  we  expect  the  scaling  of  the 
pressure  field  to  be  on  qQ  and  5*  as  outlined  earlier.  We  find  at  25  ft/sec,  6*/h  ~  0.13.  Fig¬ 
ure  16B  shows  that  separate  models  of  the  coherences  are  required  above  and  below  coh/U0 
=  3.  As  explained  at  the  end  of  Chapter  3  we  consider  that  the  low  frequency  vibration 
affects  the  boundary  layer  as  an  oscillating  angle  of  attack  which  probably  does  not  instan¬ 
taneously  exceed  2  deg. 

For  high  frequencies  coh/U0  >  3,  corresponding  to  cu5*/U0  >  5,  we  model  the  coher¬ 
ence  with  Equation  (30).  Curve  fitting  to  Figures  13  and  15  yields 


NRew  =  U0W/l; 
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Xx/6*  ~  5.5 


-l 


\J8*  ~  0.8 


Figure  14  gives  Uc  =  0.7  U0  .  The  average  of  8*  over  the  width  is  determined  by  the  down 
stream  half  of  the  width  so  that  ww/w  =  1/2,  and  £w/L  =  0.7.  Substitution  into  Equation 
(34)  with  =  2,  gives 
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where  we  have  approximated 
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Equation  (35)  is  the  high-frequency  boundary-layer  excitation  designated  in  Figure  35.  The 
estimate  agrees  well  with  the  measured  values.  Similarly,  for  the  1.5-  by  14-in.  beam,  as¬ 
suming  that  the  boundary  layer  is  nearly  similar  to  that  developed  on  the  2.75-in.  beam,  we 
have  £w  /L  =  1/2,  and  ww/w  ~  1,  which,  using  data  for  a  =  0°,  gives 
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and  these  equations  are  shown  in  Figure  3 1 .  The  theoretical  expressions  in  both  cases  ac¬ 
count  reasonably  well  for  the  observed  acceleration  levels. 

At  low  frequencies  of  coh/U0  <  3,  we  use  Equations  (29)  which  contain  \x  =  13.7  6*, 
and  Xz  =  2.9  8*.  We  also  have  Uc  —  UQ.  Further,  approximating  that 
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This  is  designated  the  low-frequency,  a  =  2°,  boundary-layer  excitation  in  Figure  35.  In  this 
case  the  observed  acceleration  spectrum  is  greater  than  that  which  could  theoretically  be  gen¬ 
erated  by  the  low-frequency,  boundary-layer  pressure.  In  the  midfrequency  range  near 
ooh/U0  ~  3,  the  measured  acceleration  spectrum  is  close  enough  to  the  theoretical  estimate 
for  response  to  this  part  of  the  pressure  field  to  be  plausible. 


RESPONSE  TO  INFLOW  TURBULENCE 

The  right  hand  side  of  Equation  (34)  can  be  interpreted  as  a  mean-square  pressure  coeffi¬ 
cient  exerted  on  the  surface  of  the  strut.  Mugridge  (1970)3  has  theoretically  determined  the 
mean  square  pressure  coefficient  on  a  rigid  strut  exposed  to  a  three-dimensional  turbulent  inflow. 
We  haveadopted  Mugridge’s3  result  to  determine  the  response  of  the  beams  to  the  water  tunnel 
turbulence.  Thus  we  obtain 

*AA<W’*0>U<>/h  ,  ,  Uo 

- : -  K  +  mrn  >  T?2  =  0L  (W)  — 

% 


where  we  have  assumed  isotropic  turbulence;  L'  is  the  macro-scale  of  the  field  and  k  =  cuw/2U0. 
We  have  assumed  L'  <  w.  This  equation  gives  the  resultant  pressure  spectra  in  terms  of  the 
inflow  turbulence  spectra  which  have  been  given  in  Figure  35.  We  have  found  L'  —  2h  and  we 
have  already  used  the  background  intensity  uB  /UQ  —  0.028. 
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For  a  strut  of  thickness  h,  we  would  expect  Equation  (36)  to  apply  for  cuh/U0  <  1  be¬ 
cause  for  higher  frequencies  the  pressure  fields  above  and  below  the  beam  would  necessarily 
be  uncorrelated.  The  result  of  Equation  (36)  is  shown  as  the  upper  dotted  line  in  Figure  36. 
The  close  agreement  with  the  measured  acceleration  levels  gives  a  substantial  indication  that 
low  frequency  acceleration  is  generated  by  the  inflow  turbulence. 

The  estimates  of  beam  acceleration  that  have  been  made  in  this  part  have  indicated  that 
inflow  turbulence  generated  the  very  low  frequency  vibration.  Boundary  layer  pressures  gen¬ 
erated  by  the  flow  on  the  beam  accounted  for  the  very  high  frequency  acceleration.  The  low 
frequency  pressures  which  are  associated  with  leading  edge  separation  may  have  accounted  for 
acceleration  only  in  a  small  frequency  region  around  cuh/U0  —  3. 


6.  CONCLUSION 

The  linear  interaction  of  a  beam  with  flow  parallel  to  its  neutral  axis  has  been  character¬ 
ized  in  terms  of  fluid  excitation  and  damping.  The  vibration  of  the  beam  at  high  frequencies 
is  controlled  by  its  own  boundary-layer  pressure,  and  the  mean-square  modal  acceleration  has 
been  found  to  be  proporational  to  both  the  correlation  area  and  intensity  of  the  pressure 
field.  At  the  frequencies  for  which  the  boundary  layer  controls  the  flow  excitation,  the 
correlation  area  is  inversely  proportional  to  the  square  of  the  frequency.  Dependence  of  the 
pressure  field  on  large-amplitude,  low-frequency  motion  manifests  itself  as  an  oscillating  angle 
of  attack  which  causes  a  somewhat  thicker  boundary  layer  than  would  otherwise  exist;  how¬ 
ever,  this  does  not  have  a  dominant  influence  on  the  level  of  beam  vibration.  The  coherence 
and  intensity  of  the  boundary-layer  pressure  field  is  controlled  by  the  flow  in  the  wedge  re¬ 
gion  of  the  trailing  edge. 

At  low  frequencies,  the  inflow  turbulence  controls  the  excitation  of  the  beam.  For 
a>h/U0  <  1,  the  response  is  given  by  the  Mugridge  3  result;  however,  for  higher  frequencies, 
existing  simple  theories  are  not  adequate  to  correctly  predict  the  response  to  inflow  turbu¬ 
lence.  It  is  reasonable  that  this  small-scale  turbulence  influences  the  response  of  the  beam 
just  as  the  relatively  large-scale  turbulence  in  the  boundary  layer  formed  on  the  beam.  Large 
eddies  in  the  boundary  layer,  such  as  those  originating  with  flow  separation  at  the  leading 
edge  of  the  strut,  could  materially  contribute  to  the  beam  vibration  at  moderate  frequencies, 
i.e.  cuh/U0  ~  3. 

The  fluid  loading  on  the  beam  is  dominated  by  inertial  and  circulatory  effects.  The 
magnitude  of  the  measured  entrained  inertia  did  not  differ  from  that  which  would  be  calcu¬ 
lated  from  simple  potential  theory;  see  Equation  (9).  Hydrodynamic  damping  was  the  domi¬ 
nant  contributor  to  the  total  loss  factors  of  most  modes  at  high  speeds.  As  predicted  by  two 
dimensional  oscillating  hydrofoil  theory,  the  hydrodynamic  damping  for  torsion  modes  was 
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double  that  for  bending  modes.  The  speed-dependent  resistance  to  vibration,  which  is  in¬ 
duced  by  the  fluid  and  which  causes  the  damping,  is  normalized  on  the  measured  entrained 
inertia  to  form  a  hydrodynamic  loss  factor.  So  defined,  this  loss  factor  is  apparently  a  func¬ 
tion  of  a  reduced  frequency,  based  on  the  width  of  the  beam  and  inflow  speed;  this  depend¬ 
ence  is,  within  the  limits  of  the  parameters  covered  in  the  experiments,  nearly  unique.  This 
observation  is  theoretically  supported,  since  application  of  oscillating-airfoil  theory  for  vanish¬ 
ing  aspect  ratio  yields  a  hydrodynamic  loss  factor  which  is  independent  of  the  aspect  ratio. 

The  vibration  patterns  measured  at  resonance  deviate  from  the  theoretically  calculated 
mode  shapes  in  vacuo.  The  deviation  is  due  to  the  nonideal  clamping  as  well  as  the  chord- 
wise  variation  in  thickness.  That  fluid  loading  does  not  affect  the  resonance  vibration  pattern 
is  indicated,  since  the  same  patterns  were  obtained  in  air  and  water. 
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Mean-Velocity  Profiles  at  various  Stations  along  the  Strut  in  the  Wind  Tunnel  at  oc  =  2  Degrees 


Figure  6  —  Turbulence  Intensities,  Measured  along  the  Rigid  Strut  at  a  =  2  Degrees 
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Autospectra  of  the  Wall  Pressure  on  the  Strut  in  the  Wind  Tunnel  at  oc  =  o  Degree 
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Figure  12  -  Normalized  Longitudinal  Cross  Spectral  Density  Magnitude 
on  the  Rigid  Strut  at  «  =  0  Degree 


Figure  13  -  Normalized  Longitudinal  Cross  Spectral  Density  Magnitude 
on  the  Rigid  Strut  at  «  =  2  Degrees 
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Figure  1 5  —  Normalized  Lateral  Cross-Spectral  Density  at  Position  C 
on  the  Strut  at  =  2  Degrees 
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Figure  16  —  Streamwise  Coherence  of  the  Pressure  Field 
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Figure  17  —  Beam  Support  and  Clamping  Arrangement 


.75- 


O CENTER  •  KNUCKLE 


OL- 


52 


OPEN  POINTS 
HALF  POINTS 
CLOSED  POINTS 


J _ L 


DRY  FRICTION 


ESTIMATED  VISCOUS 
LOSS  FACTOR  FOR 
2.75  INCH  BEAM 


1.5" 

2.75 

5" 


I  I  l _ 1_ _ I 

102 

FREQUENCY,  Hz 


J _ I — 

103 


22  —  Total  Loss  Factors  of  the  Beams  in  Water  with  U  =  0 
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Figure  24  -  Hydrodynamics  Loss  Factors  on  the  20-Inch-Long  Beams 
of  1.5-  and  2.75-Inch  Widths 
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Figure  27  -  Frequency  Spectral  Densities,  a2  /Af,  Determined  from  the  Flow-Induced 

Modal  Response  of  the  2.75-Inch  Beam 


Figure  29  -  Normalized  Acceleration  Spectra  for  the  5-  by  20-Inch  Beam 
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Figure  34  -  Autospectra  of  the  Longitudinal  Turbulence  Intensity  at  Position  F, 
x/h  =  0.2,  Measured  in  Air  and  Water 
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Figure  35  —  Spectral  Densities  of  Background  Turbulence  Intensities  in  the  Water  Tunnel 
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TABLE  3  -  COMPARISON  OF  BOUNDARY 
LAYER  PARAMETERS  IN  AIR  AND  WATER 


Air  (a  =  2°) 

Water 

w  in  inches 

22 
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50 

7.7 
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0.32 
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6/5* 
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5 

5*/h 

0.16 
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U0w 

-  NR 

v  Rew 

5.7  x  105 
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U05* 

-  nr 
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3  x  103 
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frequency  based  on  the  width  of  the  strut  and  inflow  speed.  Finally,  a  wind  tunnel  study  of  the  statistical 
properties  of  the  boundary  layer  formed  on  the  strut  is  described.  The  results  disclose  that  flow  separation 
at  the  leading  edge,  which  is  sensitive  to  angles  of  attack,  generates  a  low-frequency  pressure  field  that  is 
markedly  higher  than  that  normally  encountered  in  boundary  layers.  At  high  frequencies  the  pressure 
field  is  influenced  by  the  local  flow  parameters  normally  used  in  boundary-layer  studies. 
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